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ABSTRACT 

The present paper extends our previous theory of the stellar initial mass 

function (IMF) by including the time-dependence, and by including the impact 

Cf^ ' of magnetic field. The predicted mass spectra are similar to the time independent 

_il ■ ward smaller masses by a factor of a few. Assuming that star-forming clumps 

O ' follow Larson type relations, we obtain core mass functions in good agreement 

with the observationally derived IMF, in particular when taking into account the 

thermodynamics of the gas. The time-dependent theory directly yields an analyt- 

k>( ' ical expression for the star formation rate (SFR) at cloud scales. The SFR values 

H ■ agree well with the observational determinations of various Galactic molecular 

clouds. Furthermore, we show that the SFR does not simply depend linearly on 

density, as sometimes claimed in the literature, but depends also strongly on the 

clump mass/size, which yields the observed scatter. We stress, however, that any 

SFR theory depends, explicitly or implicitly, on very uncertain assumptions like 

clump boundaries or the mass of the most massive stars that can form in a given 

clump, making the final determinations uncertain by a factor of a few. Finally, 

we derive a fully time- dependent model for the IMF by considering a clump, or 
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a distribution of clumps accreting at a constant rate and thus whose physical 
properties evolve with time. In spite of its simplicity, this model reproduces 
reasonably well various features observed in numerical simulations of converging 
flows. Based on this general theory, we present a paradigm for star formation 
and the IMF. 

Subject headings: stars: formation — stars: mass function — ISM: clouds — 
physical processes: turbulence 



1. Introduction 

Understanding the origin of the initial mass function (IMF) and inferring the star for- 
mation rate (SFR) in galaxies are the two main challenges of star formation theory. Many 
attempts have been made along the years to try to resolve these two fundamental issues 
(see Hennebelle & Chabrier 2010 for a brief review regarding the main modern theories of 
the IMF). Within the past three years, we have developed a new analytical theory, based 
on the gravo-turbulent paradigm of star formation (e.g. MacLow & Klessen 2004), aimed 
at explaining the IMF (Hennebelle & Chabrier 2008, 2009, hereafter papers I and II). In- 
deed, in our theory of the IMF, as well as in the Padoan & Nordlund (2002) one, large-scale 
supersonic turbulence is supposed to generate small-scale overdense regions with respect to 
the surrounding background, within which gravity eventually dominates all supports and 
triggers the collapse, leading to the formation of the prestellar cores. As in Padoan & Nord- 
lund's theory, we assume that there is a direct correspondence between the prestellar core 
mass function (CMF) and the final stellar IMF. More preciselly, we assume that there is a 
good correspondance between the mass reservoir out of which the cores form and the IMF. 
A fact that seems to be supported observationally by the remarkable similarity between the 
IMF and the CMF (see e.g. Andre et al. 2010 for the most recent results concerning this 
issue) but also statistically by the rather strong correlation between the CMF and the IMF 
inferred from the analysis of numerical simulations aimed at exploring this issue (Chabrier & 
Hennebelle 2010). Recently, Hopkins (2011, 2012) derived an IMF theory in a similar spirit 
as the Hennebelle- Chabrier one, based on a different, so-called excursion set formalism, and 
has extended the results to larger scales, typical of large-scale structures in galactic disks. 
His results are found to agree fairly well with the Hennebelle-Chabrier ones when considered 
at the star forming clump scale. 

The theory developed in papers I and II, which extends to the context of star formation, 
characterized by non-linear density fluctuations, the formalism developed in cosmology for 
linear fluctuations by Press & Schechter (1974), consists in properly counting at all scales the 
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self-gravitationally bound density fluctuation^^. In the HC theory, these collapsing density 
fluctuations represent the overdense regions that isolate themselves from the surrounding 
medium and start to contract under the action of gravity at the very initial stages of star 
formation, and out of which prestellar cores, and later on individual stars (or multiple 
star systems such as binaries) will form (see Chabrier & Hennebelle 2011 for a simplified 
explanation of the HC formalism). As mentioned in papers I and II, a limitation so far of 
this theory resides in the time independent nature of its formulation. The importance of 
time-dependence for the IMF has been stressed, for instance, by Clark et al. (2007) who 
argue that, since massive (low-density) cores are expected to collapse in a time longer than 
low-mass (high-density) ones, the CMF should be flatter than the IMF, if the latter one is to 
be inherited from the former one. In the Hennebelle- Chabrier theory, however, turbulence- 
induced velocity dispersion, rather than purely thermal motions, plays a dominant role in 
setting-up the mass of the massive cores when they form, so that the turbulent Jeans mass, 
which entails a Mach dependence, instead of the thermal Jeans mass, should be used as the 
characteristic mass scale. This does not mean that turbulence is acting as a pressure or 
even a support in a static sense, as often misunderstood, but rather that on large scales, 
turbulent motions act to prevent larger amounts of material from immediate gravitational 
collapse. Eventually, for instance where the flows collide, the density is temporarily enhanced, 
yielding the subsequent collapse of these large mass reservoirs, progenitors of massive cores 
(see Chabrier & Hennebelle 2011). Under these conditions, as shown in Appendix C of 
Paper II, the characteristic free-fall timescale depends much more weakly on the star mass 
{Tff oc M^/^), making the aforementioned time problem much less severe. In the present 
paper, we develop this argument quantitatively by deriving a time-dependent theory of the 
IMF. We show that, indeed, time- dependence barely affects the slope of the IMF at large 
masses. The peak of the mass spectrum, however, is shifted by a factor of ~ 3 towards 
smaller masses, as more low-mass prestellar cores are able to form. This time-dependent 
formulation of star formation enables us to derive a star formation rate. We stress, however, 
that this SFR is valid at the scale of the clouds and not necessarily at the scale of the entire 
Galaxy (see e.g. Ostriker et al. 2010). 

A summary of the present results and comparisons with previously published SFR theo- 
ries and with some observational results have been presented in Hennebelle & Chabrier (2011, 
HCll). The present paper presents in details the whole derivation of our time-dependent 
IMF theory and of the analytical SFR and confronts the results with further recent obser- 



^As discussed in §3 and §5 of paper I, the HC theory takes into account the probabihty for overdense, 
collapsing structures to be included in larger collapsing ones, and thus properly addresses the so-called "cloud 
in cloud" problem present in the Press & Schechter formalism. 
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vational determinations. The paper is organized as follows. The extension of our formalism 
to a time-dependent derivation, and comparisons between the time- dependent and time- 
independent theories of the mass spectrum of self-gravitating fluctuations, are presented in 
details in §2, both for the isothermal and non-isothermal case. The effect of magnetic field 
is also considered in this section. In §3, we derive the star formation rate and star formation 
efficiency. We first derive the theoretical expression and then compute the SFR values for an 
ensemble of clump parameters and explore the dependence of the SFR upon various clump 
characteristic propertieqj. In §4, we make comparison with recent SFR determinations of 
molecular clouds in the Galaxy. We briefiy discuss the impact of filaments in star formation 
upon our formalism. In §5, we present the complete self-consistent time- dependent model, by 
exploring the consequences of our time-dependent theory on the mass spectrum of different 
evolving clumps. Then, we investigate the impact on the mass spectrum of a time-dependent 
clump distribution instead of a single clump. Section 6 is devoted to the conclusion and our 
paradigm for star formation and the IMF is presented. 



2. Time dependent theory 

2.1. Analytical expression of the mass spectrum 

The underlying concept on which relies the theory developed in papers I and II consists 
to identify, in a random field of density fiuctuations, the mass M^ which at scale R (physically 
speaking, R denotes the radius of the density fiuctuation) is gravitationally unstable, i.e. the 
mass contained in regions within which gravity dominates over all sources of support. To 
achieve this, the first step is to determine M^j, the mass contained in regions whose density 
exceeds a scale- dependent density threshold, log^lj, determined by the virial condition. Since 
all this gas is unstable at scale smaller than or equal to R, it is expected that it will end up in 
objects of mass smaller than or equal to the Jeans mass associated to the density threshold. 
Thus the second step is to equal Mr with the mass contained in the structures of mass equal 
to or smaller than the associated (turbulent) Jeans mass (see eqns.(28) and (29) of paper I). 

This approach implicitly assumes that each part of the fiow is initially assigned a spe- 
cific (scale dependent) Jeans mass, which collapses. While such an assumption is certainly 
reasonable in the case of cold dark matter fiuctuations within the primordial universe, it is a 
priori not the case in a turbulent fiow in which fiuctuations of scale R are replenished within 



^In our theory (see papers I and II), star forming "clumps" are identified as overdense (n > lO^-lO"^ cm"'^) 
~l-10pc-size unbound regions within large, diffuse molecular clouds, within which '^O.lpc-size gravitationally 
bound prestellar "cores" will preferentially form. 
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a few crossing times, tr ~ 2R/Vrms{R), where K-ms(-R) = {o'^y^'^ denotes the turbulent rms 
velocity at scale R and a^ the velocity fluctuation over scales smaller than R. Therefore, 
in the case of turbulent molecular clumps, it seems necessary to take such dynamical effects 
into account. This is not obvious by any means, in particular for large scale fluctuations for 
which, first, a long enough time is required and, second, gravity must not entirely freeze the 
gas motions by preventing fluid particles to escape their local Jeans masses. In this latter 
case, gravity will prevent the flow motions and the replenishment of the density fluctuations, 
which have become gravitationally unstable. 

On the other hand, if the clump is accreting at a sufficiently high rate or if it contains 
a large fraction of gas that is not dense enough to be locally gravitationally unstable (i.e. if 
the Jeans length at the gas density is comparable to or larger than the size of the clump), 
dissipation of large scale turbulence leads to turbulent compression which can continuously 
generate density fluctuations at all scales and trigger, for the densest ones, gravitational 
collapse. In this case, one must take into account the fact that, during the lifetime of the 
cloud. To, the fluctuations at scale R have been replenished a number of time equal to tq/tr. 
This implies that eq. (31) of paper I must be modified as 



Vr 



/•oo n^^R 

/ pexp((5)P/j(5)(— )rf(5= / M'Af{M')P{R,M')dM'. (1) 

is?, '^R Jo 



where 6 = log(p/p) denotes the (logarithm of) density fiuctuation, 5^ the threshold density 
at scale R, V{6) is the density distribution (PDF)tl and Vc ^ L^ denotes the clump's volume. 

Apart from the term Tq/tji, this equation is the same as the one derived in paper 
I. It uses the formulation developed by Jedamzik (1995) in the context of dark matter 
halos (see also Yano et al. 1996). The first equality for Mtot{R) stems from the fact that 
the mass contained within structures of mass M < M^ is equal to the mass of the gas 
which, smoothed at scale R, has a (logarithmic) density larger than a critical threshold 5^. 
The second expression arises from the fact that the number- density of structures of mass 
M < M^ is Af{M')P{R, M')dM'. Here, Af{M')dM' is the number-density of structures of 
mass between M' and M' + dM', while P{R, M') is the probability to find a gravitationally 
unstable structure of mass M' embedded inside a structure of gas which at scale R has a 
(logarithmic) density larger than 5^. P{R, M') is assumed to be equal to 1 (see Appendix 
D of Paper I for further justification). The time ratio tq/t/j thus simply illustrates the fact 
that the fiow fiuctuations at scale R have been rejuvenated Tq/tr times. 



■^We stress that, in our general formalism, V is not necessarily a lognormal, as emphasized in Schmidt et 
al. (2010). 
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Taking the derivative of eq. ([T]) with respect to R, we obtain for the number- density mass 
spectrum of gravitationally bound structures, M{M) = d{N/Vc)/dM: 



^'^^a - 4^(-^<S''^'"'^'^''''^' + r'^'"''^'<S'^''i''') •'' 



R 



As shown in paper I, the first term in this expression is the most important one and dominates 
over the second one except when the scale R becomes comparable to the injection scale, Li, 
which is basically the size of the cloud (see the interesting formulation of Hopkins (2012) to 
avoid this large-scale limitation of the present formalism). 

Once the expression of the critical density threshold, 6^^ = \n{p{R) / p) , is specified from 
the Virial condition, {V^^g{R)) + 3C| < Epot{R)/M (see §4.2 of Paper I), the mass spectrum 
of the self-gravitating pieces of fluid can be inferred from 

MR = CmP{R)R^ (3) 

where Cm is a geometrical factor, typically of the order of 47r/3. 

Before proceeding further, it is worth stressing what is exactly selected by our procedure. 
As seen from eq. ([T]), the integration is performed from to M^. This means that the 
smallest Jeans masses, corresponding to the densest pieces of gas, are first accounted for and 
removed from the available gas mass. Then the larger Jeans masses are progressively taken 
into account and removed. By doing so, we properly take into account the fact that there are 
small Jeans masses embedded into larger ones. This happens, in particular, when the PDF 
significantly varies locally, i.e. when turbulence is strong, which implies that in that case 
the scale must be large. Therefore, strictly speaking, we do not identify well defined bound 
cores but rather coherent mass reservoirs in the density field that isolate themselves from the 
surrounding medium under the action of gravity. For sufficiently small scales, the density 
becomes reasonably uniform within the mass reservoir to have a clear correspondence between 
this latter and a well identified "core". For large scales, which will lead eventually to the 
formation of massive stars, however, the one-to-one core-reservoir correspondence becomes 
more blurry (as the potential well associated with the mass reservoir becomes itself more 
shallow). In that case, the reservoir of mass out of which the most massive cores/stars will 
form corresponds to what is left once all the self-gravitating small-scale density fiuctuations 
embedded in the reservoir have been properly accounted for and "taken away" . In that sense, 
for the largest scales, thus the most massive cores, our formalism rejoins in some sense the so- 
called "competitive" accretion process (see e.g. Smith et al. 2009). Therefore, in principle, 
the turbulent fragmentation process is properly described in our formalism, apart from the 
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fact that, for sake of simplicity, the conditional probability P{R,M) is taken to be equal to 
1. What has not been taken into account so far, however, is the gravitational fragmentation 
that occurs during the collapse. Although such fragmentation may occur, its importance 
is likely to remain limited because of the impending roles of magnetic field (Machida et al. 
2005, 2008, Hennebelle & Teyssier 2008, Price & Bate 2008, Hennebelle et al. 2011) and 
radiation (Krumholz et al. 2007, Bate 2009, Commercon et al. 2010) particularly when both 
are present (Commercon et al. 2011), an issue which still needs to be properly quantified. 
Observations of massive cores indeed suggest that fragmentation is rather limited, most of 
the mass of the core ending up in one or just a few smaller cores (Bontemps et al. 2010, 
Longmore et al. 2011, Bressert et al. 2010, Palau et al. 2013). 



2.2. Influence of the time-dependence on the CMF/IMF: isothermal case 

In this section, we examine the impact of time-dependence on the CMF in the simple 
isothermal case. The barotropic case will be examined in §2.31 We first discuss the density 
PDF, then the crossing time. Finally, we derive the analytical expression and compare the 
results with the time-independent ones. 



2.2.1. Density probability function 

Numerical simulations (Vazquez-Semadeni 1994, Padoan et al. 1997, Passot & Vazquez- 
Semadeni 1998, Kritsuk et al. 2007, Federrath et al. 2008, Schmidt et al. 2009) have 
revealed that the turbulence- induced density distribution, V{6), is reasonably well described 
by a lognormal distribution 

P(5) = — ==exp' ^ ' 



S = Hp/p), 5 = -al/2, al = ln(l + b'M'). (4) 

In this expression, Ai is the characteristic cloud Mach number and b a non-dimensional 
coefficient that depends on the turbulence forcing (see Federrath et al. 2010). It typically 
varies from 0.25 when the forcing is purely solenoidal to almost 1 when the forcing is applied 
only on compressible modes. Such a 3D lognormal shape for density fluctuations has received 
observational support from its 2D observed projection, namely the power spectrum column 
density of molecular clouds, measured from dust extinction maps (Kainulainen et al. 2009, 
Brunt et al. 2010). The scale-dependence of the variance of the distribution reads, in 3D 



(see Paper I) 



a\R) = I SikfA-KkHk = al\l-(^] 1 , (5) 



2-k/Lc \ V Lr 



where L^. is the cloud's size, 5(/c)^ oc A;^" is the power spectrum of log(p), of 3D index n' . This 
latter is found in isothermal, shock-dominated hydrodynamical and MHD simulations to be 
very similar to the index n of the velocity power spectrum, with a typical value n' ^ n ^ 3.8, 
between 11/3 (Kolmogoroff limit) and 4 (Burgers limit) (see e.g. Kritsuk 2007, Federrath et 
al. 2008, Schmidt et al. 2009). Then (see paper I) 
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2.2.2. The crossing time, tr 

The crossing time at scale R is given by tr = Tct{R) = 2actR/Vct, where Vd is the 
relevant velocity and a^ a dimensionless coefficient of the order of a few. At large scales. 
Vet is typically equal to the one-dimensional velocity dispersion l^|^ = Kms/\/3, where 
l^nig = P^^ designates the 3D velocity dispersion all along this paper, while at small scales, 
below the sonic length, V^ ~ Cg- This crossing time is the typical time that is necessary 
for the density field to be significantly modified at scale R, implying that a new set of 
fiuctuations, statistically independent of the former one, has been processed. 

As mentioned earlier and discussed in papers I and II, in the present context, we select 
the pieces of gas which are self-gravitating, i.e. such that their internal gravitational, kinetic 
and thermal energies obey the condition: — -Egrav > 2i?kin + 3Pt/i- At large scales, this implies 
that agGM/R > V^^s) where ag is a dimensionless coefficient, equal to 3/5 for a uniform 
density fiuctuation. We thus get 






^R = -THd- = 2arf J^^:^r//(i?), (7) 



rms 



where Tff{R) = J z2G{r) ^^ ^^^ free-fall time of the density fiuctuation of density p{R). A 
similar expression is obtained below the sonic length. In the following we thus define the 
crossing time of a collapsing density fiuctuation of scale R as 

TR = (PtTffiR), (8) 

/ 24 
with 0i = 2aetW^^3. (9) 
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Defining r°y as the free-fall time at the clump's mean density, r°j = J ^p — l-0'<'(2^)~^^^( io3cm-a )~^^^ 
Myr (where // = 2.33 is the mean molecular weight for a cosmic H2/He composition), we get 

-^ = <Pt\-- (10) 

Note that the choice of r9r is not consequential at this stage as it simply modifies the value 
of A/" uniformly without affecting its shape. Combining eqs. ([2]), (|6]), and ([9]) and dropping 
the second term in eq. ([2]), as mentioned earlier, and assuming that r^j ~ tq, we obtain (see 
paper I) 

(~ X 1/2 ~ 

i?3 ) dM% dR V2^ "^ V 2^2 2 8 ; ' 
where M]j is the critical mass at scale R (see eq.(j3])) while 

^« -!.>(§). (12) 

where R = R/X^j, M = M/Mj and Mj, Aj denote the usual mean thermal Jeans mass and 
Jeans length, respectively: 

where a j is a dimensionless geometrical factor of order unity. Taking the standard definition 
of the Jeans mass as the mass enclosed in a sphere of diameter equal to the Jeans length, 
one gets aj = 7r^/^/6. We now need to know M]^, and its derivative dM^/dR, as a function 
of the scale, R. It is determined by the physical processes at play in the cloud, as examined 
in the next sections. 



2.2.3. Analytical expression 

We now need to specify the density threshold, 5^, which is determined from the virial 
theorem. As shown in paper I (§4.3), the condition for collapse, which simply reads 5 > 5'j^ 
or equivalently M > M"^ = Mj{R), yields after calculations 

T/2 1 3/2 

(c.(p))^ + (V 



M > aj- 



v/G'3pexp(5) 



(15) 
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We assume that the effective dispersion velocity obeys the Larson's (1981) relationship: 



(KL) = K,'x 



Ipc 



M = M{R) = ^ '"^^ . (16) 

As mentioned in paper I, the Mach number Ai represents "effective" values which include 
both the hydrodynamical and magnetic contributions, i.e. Kms = {(K-ms)hydro + ^4/2}^^^ 
where Va = B /{Aixpy/'^ denotes the Alfven velocity. The coefficient r] is related to n, the 
index of the velocity power spectrum, by the relation (see eq.(24) of Paper I): 

72 — 3 , , 

^=^-- (17) 

As mentioned earlier, 3D numerical simulations of compressible turbulence (e.g. Kritsuk et 
al. 2007) suggest a value n ~ 3.8, yielding rj ~ 0.4, as indeed found in observations. 

With eq.(j3]), eq.( IT5|) implies (see Paper I) 

M>M% = af{^R+^^, (18) 

After normalisation, eq. (I18p becomes (see §5.4 Paper I) 

M^ = M/M^j = R{1+ MlR^"^), (19) 

where A^* is given by: 



and thus illustrates the impact of turbulence induced velocity fluctuations at the Jeans scale 
(see paper I). 

Equations (ITTj) and (IT9l) finally yield for the mass spectrum of gravitationally bound 
prestellar cores: 

M(M) = '-K i ^^"-'"■^T X M] " X -P<--'/^'. (21) 

where A/q = p/Mj. Equation (|2T|) is the time-dependent generalization of eq. (44) of paper 
I. The time dependence appears explicitly through the factor 1/0^ but also through the 
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modification of tlie exponent, —1 — ^\n{M/R^), instead of —3/2 — ^ln(M/R^), wfiicfi 
arises from the time correction, proportional to y/p oc [M / B?)^/"^ . As expected from the 
discussion in the introduction, the time dependence thus affects the slope of the CMF. We 
quantify this effect in the next sections. 



2.2.J1.. Results and comparison with the time independent model 

Before comparing the time-dependent and time-independent distributions, we determine 
the position of the peak of the CMF. As discussed in paper I (§7.1.4), the peak occurs in 
the thermally dominated regime, which corresponds to A^^ = 0. The derivative of eq. fl^T]) 
with respect to the mass yields in that case 

M,».^^^exp(-.^)^i^^, (22) 

whereas in the case of a time independent distribution we have (eqn.(46) Paper I) 

Mpeak = exp(--a2) = (^^^2_^2)3/4 - (23) 

This implies that, for a given Mach number, the peak of the distribution is shifted toward 
smaller masses when time dependence is taken into account. This is intuitively expected 
from the fact that more small-scale objects, with short free-fall times, are continuously 
produced during the clump evolution in a time-dependent collapse. It is thus important to 
stress that in the (both time-dependent and time-independent) HC theory of gravo-turbulent 
fragmentation, the peak (i. e. the characteristic) mass of the IMF depends not only on the 
clump's mean thermal Jeans mass, hut also strongly (cc A4.~''^) on the characteristic Mach 
number. 

Another important quantity is the exponent of the high-mass tail of the distribution. 
Indeed, in the limit M ~ Al^i?^^"*"^, which corresponds to the turbulence dominated regime 
(see eq. dH])), eq. f l?T]) tends to a power law plus a lognormal, such as A/" oc M"*^^"*"^) with 

whereas in the time-independent case, this coefficient reads (see eq. (43) of paper II) 

As expected, the distribution becomes slightly steeper when the time dependence is included 
in the derivation of the CMF/IMF. This is due to the fact that some massive clumps fragment 
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Fig. 1. — Isothermal case. Core mass function, dN/d\ogM, for A^ = 6, A^* = v2 and 
mean clump density n = 3 x 10^ cm^^ (top panel) and n = 5 x 10^ cm"'^ (bottom panel). 
The solid line corresponds to the time-dependent model while the dotted line represents the 
time independent one. The dashed line is the Chabrier IMF shifted upward in mass by a 
factor 3 (see text). 
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into smaller pieces during the cloud collapse. Typical values for x range from 1.1 to 1.5, 
depending on r], a and A^^,, bracketing the Salpeter value, x =1.35. 

Figure d] displays the core mass function, dN/dlogM for r] = 0.45, A^ = 6, Ai^: = \/2 
and two typical clump densities, n = 5000 cm~^ and n = 3 x 10'^ cm^^. Solid lines display 
the time dependent results, dotted lines the time independent ones. For reference, dashed 
lines represent the Chabrier system IMF (Chabrier 2003), shifted upward in mass along the 
X-axis by a factor 3 to account for the observed shift between the CMF and the IMF. In 
the rest of the paper, we will refer to this IMF as the "shifted Chabrier IMF" (SCIMF). As 
mentioned above, the time-dependent distribution peaks at lower masses and has a slightly 
steeper high-mass slope than the time-independent one. In spite of these differences, both 
distributions match well the SCIMF above about the mean Jeans mass. At low masses, 
both distributions appear to be too narrow. As already discussed in papers I and II, this 
is essentially due to the isothermal approximation for the equation of state of the gas (see 
Paper II and below), but also to the assumed exact correspondence between the initial core 
mass and the final star mass. Indeed, as shown in Chabrier & Hennebelle (2010), taking 
into account some (expected) dispersion between the mass of the parent core and the one 
of the final star, i.e. between the CMF and the IMF, naturally leads to a broadening of the 
latter one compared with the former one in the low-mass regime. At last, gravity may also 
broaden the PDF. 



2.3. Barotropic equation of state and magnetic field 

Here we compare the time dependent and independent models when a barotropic equa- 
tion of state is used. We also introduce the magnetic field in our model and investigate its 
impact on the mass spectrum. 



2.3.1. Formalism 

Again, assuming that the pieces of fluid which collapse are gravitationally dominated, 
eq. (iT5l) becomes 



M> aj- 



3/2 

{Cs{p)f + (KL/3) + (KV6) 



v/G3pexp(5) 



(26) 
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where Va is the mean Alfven speed. Based on observational and numerical results, we assume 
that the magnetic intensity correlates with the gas density as 



P 



lb 



(27) 



Typical values for Bq and p are 10-20 /iG and 10^ cm""^, respectively (Crutcher 1999). The 
coefficient 7^ seems to depend on density. In the diffuse gas, 7;, is typically smaller than 0.5 
and around 0.1-0.3 (Troland & Heiles 1986, Hennebelle et al. 2008) while for higher densities 
when the gas is self-gravitating, it is of the order of 0.5 (although see Crutcher et al. 2010). 



After normalisation, eq. ( 126|) becomes 



M^ = M/M^j = R 






+ MlR^'^ + iV:) 
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(2^ 



where 



V* 



Bn 



VgV^Cs 



(29) 



As emphasized in paper II, the thermodynamics of the gas has a drastic impact on the 
low-mass end of the mass function. As in paper II, we use the barotropic eos suggested 
by Larson (1985, see also Glover & MacLow 2007) from observations of molecular clouds, 
namely: 



T oc p^i-^ p < p'^"*. 



T (X p 



,72-1 



p> p 



crit 



(30) 



with 7i ~ 0.7, 72 — 1 — 1.1, and p'^"'^ corresponds to the density above which dust grains 
become thermally coupled with the gas and thus dust cooling becomes the dominant cooling 
mechanism instead of line cooling. Observations suggest that p 



crit 



10 -^^ gem ^, 



I.e. 



n 



crit 



2.5 X 10^ cm ^. We keep the same eos than the one we used in paper I, inspired from 



Larson (1985) and write 



C 



- / \ m 

Ucin^_r-A + 
\ p J 


/ \ m- 

\ p J . 


l/m 


cir 


. N (7i-l)m 

- + 
\PJ 


[ciJ [pj 



- 1/m 



(31) 
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where m is a real number of order unity. In the following, m will be equal to 3 (see Fig. 4 
of paper II). Typically, 71 ^ 0.7 while 72 ^ 1 — 1.1. The critical density, p"'**, at which the 
transition between the two regimes is occurring, is expected to be about 10^^^ gcm~^ and 
we define 

In order to get the mass spectrum from eq. (fTTI) . we need to know M^ and dM'^/dR 
as a function of R. While for the first, it is not possible to get an explicit relation from 
eq. (125]) . the second can be obtained by deriving eq. ( 125]) with respect to R. The analytical 
expressions that we obtain are identical to eq. (38) of paper II except for the extra terms 
related to the magnetic field. We write them explicitly for completeness. 



dM% B , , 



i?3 dp 

C = 1-R^^—^, (34) 

dp 

D = (p<^^-^)" + (ircrit)'"p^^'"'^™)'^" + (K*)V"'"'- 



2.3.2. Comparison between time dependent and independent mass spectra in the barotropic 

case 

When a barotropic eos is considered, the mass spectrum depends on the clump's mean 
density not only as a normalisation factor (as for the isothermal eos) but also through the 
eos itself. As a fiducial value, we adopt a critical density (eq. f l32|) ) n'^"^ = 2 x 10^ cm~^, i.e. 
pcrtt ^ lo^is gcm~^, as inferred from observations. The density and velocity dispersions are 
taken to follow Larson (1981) (see also Falgarone et al. 2004) type relations 

/ D \ -0.7 
n = (do X 10^ cm-3) , (35) 

Vipcy 

^rms = (Mox0.8kms-^) (— j . (36) 
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Fig. 2. — Barotropic and magnetized mass spectra for clumps of size Re = 0.5, Re = 2, 
Re = 5 and Re = 20 pc. The clump density and Macli numbers are given by Larson-type 
relations (eq. ( I35l) ) with uq = 1 and do = 5,4,3,2 from top to bottom. Left panel: the 
solid line corresponds to the time-dependent CMF while the dotted line represents the time 
independent one. The dashed line is the Chabrier IMF shifted by a factor 3 in mass (see 
text). Right panel: the solid line is identical to the left panel. The dotted line corresponds 
to V^ = 1, 7fc = 0.1, the dashed line to V^ = 3, ■jb 
lb = 0.3. 



0.1 and the dot-dashed one to V^° = 1, 
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According to eq. flT6l) . Vq = uq x 0.8 km s~^. Figure [2] (left panel) shows the mass spectrum 
for four clumps of size Re =0.5, 2, 5 and 20 pc. The value of Uq is taken to be equal to 
1 in all cases while the top, middle and bottom panels correspond to do =5, 4, 3 and 2, 
respectively. The four panels thus correspond, from top to bottom, to typical clump masses 
and mean densities, Mc=183, 3564, 21995 and 355607 Mq, n = 8120,2460, 970 and 240§. 
As in Fig. [H solid hues display the time dependent case, dotted lines the time independent 
one while the dashed hne corresponds to the SCIMF. In the four cases, except possibly for 
the smallest clump (as expected, see §6.2.2 of paper II and §3. 2. 2b below), the agreement 
with the SCIMF is very good. Importantly enough, the clump densities found to yield a 
good agreement with the SCIMF are smaller by a factor of about 2-5 than what has been 
found in the time-independent case for similar agreement (Fig. 8 of paper II) and thus agree 
fairly well with the "standard" Larson density normalization values, c?o — 3. This is a direct 
consequence of the peak position being shifter toward smaller masses, as discussed earlier. 
Another noticeable improvement is the width of the distribution, which is larger than for the 
isothermal cases presented in Fig. [H As discussed in paper II, this improvement stems from 
the larger compressibility of the gas for a softer than isothermal eos, promoting small-scale 
collapsing structures, and from the strong dependence of the peak position on 7 (see fig. 3 
of paper II), which becomes even more acute in the time dependent CMF, whose peak is 
shifted toward smaller masses compared with the time-independent solution. 



2.3.3. Influence of the magnetic field on the core mass function 
To study the influence of the magnetic field on the mass spectrum, we define 

where -Bg'^ is a reference magnetic intensity equal to lO/xG, n is the cloud mean particle 
density and p = rripfi. We consider four cases, hydrodynamical (which serves as a reference), 
■jb = 0.3 with V^ = 1, and •jb = 0.1 with V^ = 1 and V^ = 3. We investigate primarily values 
of 76 smaller than 0.5 because our analysis consists in analyzing the density fluctuations 
generated by supersonic turbulence and to identify the ones which are self-gravitating and 
that subsequently will be amplified by gravity. Therefore, at the stage where the analysis is 
performed, the correlation between magnetic intensity and density is the result of turbulent 
processes rather than gravity, implying B ex p0i~0-3. Note that, as emphasized in Hennebelle 
& Chabrier (2008), the case jb = 0.5 is formally equivalent to the isothermal case and all 
the magnetic results can be obtained by a simple renormalisation of the rms velocity and 
Mach number (see eqn.(16)). 
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The impact of the magnetic field on the mass spectrum can be inferred from a com- 
parison between the various magnetized cases and the hydro dynamical one (Figure [21 right 
panel). When V^" = 1 and % = 0.3 (dot dashed line) the mass spectrum is essentially shifted 
toward larger masses with respect to the purely hydrodynamical case. This stems from the 
magnetic support, as expressed by eq. (ITSll . which adds up to the thermal pressure. Indeed, 
for this value of ■jb, the magnetic pressure density dependence, Pmag oc p°'^, is very similar 
to the thermal pressure density dependence, in particular because the Larson eos has an 
effective exponent at low density which is equal to 0.7. 

When 7b = 0.1, the impact of the magnetic field is different. For V^ = 1 (dotted line), 
the departure from the hydrodynamical case (solid line) is only marginal. This is because the 
magnetic support increases very slowly with density; thus at high density (which corresponds 
to low mass cores), the magnetic support is less important than the thermal one. At high 
masses, on the other hand, the support is dominated by turbulence. For V^ = 3 (dashed 
line), the distribution tends to be shallower at high masses. As discussed in paper II, the 
exponent of the mass distribution at high masses, J\f{M) oc M~^^~^^\ when only pressure 
terms are included (A^^ = 0) is given by (see eq. 32 of paper II) 1 + x = (9 — 67)/(4 — 37), 
so that if 7 < 0.2 the slope becomes shallower than the Salpeter value. 

Overall, for a typical magnetic field intensity in the diffuse (low density) gas of the order 
of S ~ 10 /iG, and a density dependence such that 7^ < 0.3, the mass spectrum is not very 
different from the hydrodynamical case. For 7^ = 0.3, the CMF still resembles the SCIMF but 
is shifted toward larger masses, a direct consequence of the magnetic support, which shifts the 
characteristic Jeans mass to larger values. When the magnetic field intensity is significantly 
larger than the aforementioned value (e.g. 30 /xG, like in the considered exemple), however, 
the mass spectrum becomes too shallow at large masses and starts departing significantly 
from a Chabrier/Salpeter IMF in this regime. Therefore, one of the predictions of our theory 
is that diffuse magnetized environments with magnetic field intensities largely exceeding 
about 10 fiG should have high-mass IMF tails shallower than the Salpeter value. 



3. Star formation rate 

In this section, we derive the SFR from our time-dependent analytical theory of star 
formation described in the previous section, analyze its dependence upon various parameters 
and clump properties and compare the results with recent observational determinations. 
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3.1. Theoretical considerations 

Following Krumholz & McKee (2005, KMK), we define the dimensionless star formation 
rate per free-fall time, SFRfj, as the fraction of clump mass converted into stars per clump 
mean free-fall time, i.e. the free-fall time defined at the clump mean density p: 

SFRff = ^r'ff. (38) 



M, 



c 



where M* denotes the total star formation rate arising from a clump of mass Mc and volume 
Vc — Rl- The star formation efficiency (SFE), is defined as the global mass fraction of a 
clump converted into stars during the lifetime Tq of the clump, which can last a few free-fall 
times (e.g. Murray 2011): 

SFE = ^ = SFRfj X (-^). (39) 

According to these definitions, SFRff is thus given by the integral of the mass spectrum 
specified by eq. ([2]), as 



SFR 



if 



Mc 



MM{M)dM 



P 

^^^"* dR ( ddRT^ff ,, ^^ ,, ^ f°° d ^° 



Tf 



dM \ dR tr Js^ dR tr J 

In this expression e is the (supposedly uniform) efficiency with which the mass within the 
collapsing prestellar cores is converted into starcl. Indeed, during the collapse, a substantial 
fraction of the mass initially within the core is blown away by outfiows and jets. Calculations 
(e.g. Matzner & McKee 2000, Ciardi & Hennebelle 2010), as well as observations, suggest 
that 

e~ 0.3 -0.5. (41) 

In eq. (140!) . Mcut corresponds to the mass of the largest star that can possibly form in the 
cloud. Equivalently, according to the mass-scale relation given by eq. (1191) or (128 p . one can 
define i/cut = R/Rc, which denotes the largest size fiuctuations that can turn unstable in 
the cloud. The value of Mcut is crucial in setting SFRff. Indeed, as emphasized in paper 



''Note the typo in eq.(7) of Hennebelle & Chabrier (2011) where it should simply be dM and not ^^ in 
the integrand. 
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I (§5.1.3), the integral Vc J^ MAf{M)dM is equal to the mass of the system itself because 
integrating up to infinity implies that all pieces of gas, including the very diffuse ones, are 
Jeans unstable. Therefore, restricting the integration to a finite value implies that any piece 
of fluid which is not dense enough, i.e. whose density corresponds to a Jeans length larger 
than a significant fraction of the cloud's size, is excluded from the mass spectrum of collapsing 
structures. The immediate consequence of such a truncation of the integral is obviously to 
reduce the SFR. 

If Re is the clump's radius, more generally its characteristic scale, the question is thus 
to determine up to which fraction of this scale or up to which fraction of the clump mass, the 
integration should be performed ? A similar question concerns the behaviour of the density 
power spectrum whose scale-dependence is characterized by cr^(-R) (eq. E]). While in the 
inertial domain, it is well established that the power spectrum of log(p) exhibits a power law 
behaviour (Beresnyak et al. 2005, Federrath et al. 2008, Audit & Hennebelle 2010), this may 
not be the case when the scale R is approaching the clump size Re. This question is directly 
related to the very definition of the clump itself and how it connects to the surrounding 
medium. Unfortunately, these questions are far from being settled. Finally, as discussed in 
paper I, the second term of the right hand side in eq. (I40j) is important only when i/cut -^ 1- 
Therefore, given all the uncertainties when approaching this limit and because the influence 
of this second term remains limited, we elected to drop it in the rest of the calculations. In 
that case, after proper normalisation, eqs. ( ITOj) . ( fT2l) . (fTTl) and ( l40l) yield 



SFR^ff = - \^^\^^ S , exp -^ ]dM, 42 

" Jo \R^ I dM-^ dR V2^ "^ V 2a2 8 ' ^ ' 




with 



SFRff = (-^) X SFR]f. (43) 



For the case of an isothermal eos, for which the mass-size relationship is given by eq. flT^ . 
this yields 



1 l„^ A/ 



R^) V2^a 

(44) 

Conversely, the SFR can also be derived from eq. ( HOj) (still dropping the second term) as: 
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(45) 



We stress that i/izs zs possible only because the second term in eq. ^B^j has been dropped. 
The value of ]}cut in eq- p5|) can be derived from the collapse condition derived in our theory 
(see eq.(29) of paper I), namely p = e^^ > e^R, or similarly from the mass-size relations given 
by eq. ( 1T9|) or (|28|) . For the isothermal case, this simply yields 



Pcut 



R-^,a + MlR'Z), (46) 



where Rcut = iVcutRc/^^j)- As clearly expressed by this equation, physically, pcut is the min- 
imum density for which a perturbation, whose associated total (thermal -|- turbulent) Jeans 
length is equal to a maximum fraction i/cut of the clump's size, can lead to a gravitational 
instability. Indeed, eq. (l46p can be rewritten: 



Mcut = ^-^^ (47) 

It is important to stress that this procedure differs not only quantitatively but also qual- 
itatively from the ones defined in KMK and PN. In these author formalisms, the critical 
density corresponds to a new, arbitrary scale, respectively a sonic scale or a shock scale, 
characteristic of a process supposed to be necessary for the collapse. In contrast, in our 
formalism, the only relevant scale is the size of the cloud itself. So, even though this can 
be expressed as a density in eq. (H5|) . the proper way to look at it, as clearly expressed by 
eq. (H7|) . is in terms of the maximum size of the fluctuations, for a given cloud's size and 
density, which are able to grow and lead to gravitational instability. Moreover, in the KMK 
and PN formalisms, the critical density is proportional to A^^, implying that only very dense 
structures will lead to star formation. In contrast, in our formalism, any structure of any 
density can collapse, if its gravitational energy dominates over all sources of support, as long 
as the associated perturbation can grow and become unstable. In practice, our results do not 
depend significantly on the value of ycut, except when y^ut ^ and for low Mach numbers 
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(see next section). Indeed, Hcut only affects the limit of the integration, which corresponds 
to a regime where the (log normal) PDF is small. 

As defined above, SFR^rr represents the SFR calculated from our theory for a core-to- 
star mass conversion efficiency e = 1 during the collapse, and for a characteristic time within 
which new mass reservoirs M of scale R become gravitationally unstable, i.e. new cores are 
produced, equal to one free-fall time at the core's density, i.e. tr = Tff{R) {(pt = 1). In that 
sense, SFR^j-f represents the " core formation rate" per free-fall time. As already mentioned, 
observations and simulations point to e ~ 0.3 — 0.5 while our estimate of Tpt in §2.2.21 leads 
to (^t ~ 3. Consequently, we typically have 

SFRff 



SFR'ff 



0.1-0.2. (48) 



3.2. Results. Dependence upon clump properties 

3.2.1. Models and assumptions 

A detailed comparison with the theoretical SFR's derived by Krumholz & McKee (2005, 
KMK) and Padoan & Nordlund (2011, PN) has been presented in Hennebelle & Chabrier 
(2011, HCll) and thus will not be repeated here. As mentioned in this paper, the key 
differences between these two theories and the present one are twofold. First, both the 
KMK and PN theories assume that the SFR is determined by one single typical free-fall 
time, defined as the free-fall time at the cloud mean density, r9r. Given the very clumpy 
nature of molecular clouds and the enormous range of density fluctuations present within 
these entities, defining one single mean free-fall time for star formation seems to be hardly 
justified. This point becomes particularly accute given the fact that there is no evidence that 
most of the regions within molecular clouds are collapsing, and if so that they are collapsing 
at r9r (see e.g. Kennicutt & Evans 2012). The second essential difference resides in the 
fact that both KMK and PN assume an ad-hoc critical density pcrit-, for star formation to 
occur (see HCll). Their SFR is thus simply obtained by estimating the gas fraction with 
density larger than Pcrit- In contrast, in our theory, the free-fall time density dependence 
of any collapsing structure of density pr is properly accounted for as the crossing time tr, 

— 1/2 

proportional to p^ , consistently varies with M and R. Therefore, a dense core of density 
p collapses in a density- dependent free-fall time Tff{p), which can differ significantly from 
the clump's mean free-fall time. In our theory, there is no particular scale or critical density, 
as we sum up over all gravitationally unstable overdense regions, whatever their scale or 
density (see eqns. ([I])- ([2])). This is indeed expected if the density fluctuations are induced 
by turbulence, which is by essence a multi-scale process, and by the fact that any piece of 
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fluid can collapse if its (turbulent+thermal) kinetic energy is dominated by its gravitational 
energy. 

Comparisons between eqns.fH^. f HS]) and the SFR's of KMK and PN have been pre- 
sented in HCll and will not be repeated here. It was shown in this paper that the KMK 
SFR differs by more than one order of magnitude from the other ones, underestimating the 
observed SFR by a similar amount, as already noticed by Heiderman et al. (2010). The 
PN SFR and the ones given by eqns. fl42l) and fl45p agree within a factor ~2-3, provided the 
(e/0t) factor is properly included in the PN relation (see HCll). The main reason why 
the present determinations yield a larger SFR than KMK|3 and PN is that, when properly 
accounting for the density-dependence of the collapsing structure free-fall times, denser fluc- 
tuations, in particular the ones denser than a given arbitrarily defined density Pcrit, collapse 
faster than lower density ones, leading to a larger formation of low-mass cores. Although 
qualitatively similar to the PN results, however, ours still differ from these latter (see Fig. 1 
of HCll). Indeed, while SFR'jr tends to decrease steadily with increasing virial parameter, 
Ovir = 2-E'kin/-£'grav, in the PN theory, it exhibits a more flatfish (nearly constant) behaviour 
before decreasing more steeply at high a^i^ in our case, as will be examined in details below. 
The reason stems from the gas more diffuse than pcru not being taken into account in the PN 
model while being accounted for, if indeed collapsing, in our model. As mentioned above, 
in our theory, if a piece of fluid of even very low density (e.g. smaller than pcrit in PN) has 
a size significantly larger than the (thermal or turbulent) Jeans length at its density, it is 
subject to collapse and thus must be taken into account. This point will be crucial when 
comparing to SFR observational determinations at low density (§!]). As mentioned earlier, 
this is a fundamental difference between our and both KMK and PN theories. 



3.2.2. Star formation rate: results 

We now examine the dependence of our theoretical SFR upon various clump character- 
istic parameters and physical properties. 

(a) Dependence upon the clump's large-scale cut-off, size and density 

Figure [3] portrays SFR^rr, for the non-isothermal case, for b = 0.5 (left column) b = 0.25 
(right column) in eq. (|1]), as a function of the cut-off parameter, i/cut = R/Rc, for various 



^besides the fact that KMK further assume that the critical free-fall time is equal to the free-fall time of 
the cloud, which implies that Pcrit = 1, a very consequential issue, as shown in HCll. 
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Fig. 3. — Value of SFR^rr for various clump sizes and densities, namely (from left to right) 
do = 4.5,3,1.5,0.7,0.1 in the non-isothermal case. Left column: b = 0.5, right column: 
b = 0.25. 
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clump sizes and mean densities (this means that the parameter do, which determines the 
density normahzation at 1 pc in the Larson relation (eqn. f l35|) ). is constant along each line 
but varies from one line to an other). As mentioned earlier and expected from the truncation 
of the integral in eq. f H2]) . SFR'jj: is basically independent of ycut for |/cui> 0.3 but quickly 
drops as this parameter approaches very small values. As also expected, the strong decrease 
of SFR^jrj occurs at higher i/cut values for smaller, less massive clumps and the truncation 
becomes more drastic as density decreases. Indeed, as already noticed in papers I and II 
(see Fig. 8 of paper II), the core mass spectrum is truncated both at high and small masses 
compared with the observationally determined CMF/IMF for small ( < 0.5 pc) clumps, with 
a drastically reduced CMF for low-density clumps. This stems from two reasons: (i) the 
scale-dependence, based on Larson's relations, of both the global and local Mach numbers, 
Ai and Ai^,, which enter our formalism (see eqns. flTB|l and (120]) ) and are responsible for 
generating and stabilizing the initial density fluctuations leading eventually to prestellar 
cores. The smaller the clump the smaller these values and thus the narrower the core mass 
spectrum (see paper I); (ii) the fact that for small clumps or low-density clumps, the Jeans 
scale becomes comparable to or larger than the size of the clump itself (see eq. (fl^I) , inhibiting 
gravitational collapse. 

According to these results, star formation is thus predicted to occur dominantly in the largest, 
most massive clumps and/or in the densest (parts of) clumps, a conclusion indeed supported 
by observations. As seen in the figure, for densities n>2000 cm~'^ and for ?/cMt> 0.1-0.3, 
depending on the clump's size, the SFR depends only weakly on ycut and reaches values 
in the range SFR^rr ~ 0.3-3.0, depending on the size of the clump, for the b = 0.5 case. 
The fact that SFR^rr, the SFR per free-fall time, can be larger than 1 is due to its density 
dependence. Indeed, SFR^j-^ is defined with respect to rj?j, the free-fall time at the cloud's 
mean density, but fluctuations of size R whose density is larger than p collapse in a shorter 
time. As already discussed, the value of ycut beyond which no star can form is ill-determined 
and may depend upon cloud parameters. In the rest of our calculations, we will assume 
ycut — 0.1 — 0.3 (see previous section), which means that only perturbations whose size is at 
most of the order of about one tenth to one third of the cloud's size are relevant to produce 
gravitationally bound prestellar cores. As just mentioned, for larger values of ycut, the SFR 
remains basically unchanged. 

As seen from Fig. [3l low-density clumps (n<1500cm^^) smaller than about ~ 2 pc 
yield quite small or even negligible values of SFR^^j-, except for values of ycut approaching 
unity, a rather unlikely possibility. Phrased differently, at low density very large clumps, but 
only very large clumps still contribute appreciably to star formation. As mentioned in the 
previous section, this is one of the direct consequences of our theory and can not be the case 
in a theory arbitrarily defining a density threshold for star formation. A point of importance 
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when we will compare with observational determinations (§!]). As just mentioned, star 
formation should thus occur dominantly in dense enough (regions of) clumps or in very 
large and massive clumps and to increase with clump density and clump mass/size. It is 
worth stressing that this strong dependence upon density or mass for efficient star formation 
naturally arises in our formalism and does not stem from an ad-hoc threshold condition for 
star formation. 

Finally, the figure shows the dependence of the SFR upon the coefficient b entering 
eq.(jl]), which relates the rms Mach number to the width of the density PDF, o"o (see eq.(jl])). 
Although the results are qualitatively similar, globally, SFR^rr is smaller by a factor of a few 
for large values of i/cut and drops drastically below larger values of ycut for b = 0.25 than for 
b = 0.5. Indeed, as shown in details in §3.1 of paper II, the mass spectrum peaks at larger 
masses and the density PDF gets narrower for smaller values of b, because the gas is less 
compressible. Since the respective contributions of the turbulence forcing mechanisms (and 
thus the value of b) in molecular clouds are rather ill-determined and may vary from cloud to 
cloud, depending for instance on the presence of expanding HII regions around massive stars 
or supernovae explosions, this implies significant possible variations of the SFR, depending 
on the environment. 

(b) Dependence upon the virial parameter 

The dependence of the SFR upon the characteristic virial parameter a^r of a clump of radius 
Re = Lc/2 and mean density p, a^r = '^Ek/Eg = {5 / 3n)V^^^/ {nG pL^) , which measures the 
ratio of turbulence over gravitational energy within the clump, has been examined in HCll 
in the isothermal case and will be examined in the next subsection in the non-isothermal 
one. It is illustrated in Fig. |U which displays the values of SFR^rr for various clump 
sizes/masses, for three typical Mach numbers, A^=16 (top), 9 (middle) and 4 (bottom). As 
a general trend, as discussed in HCll, increasing a^r leads to a decrease of the SFR, with 
an abrupt reduction above some typical value of ayir, which depends on the Mach number. 
This stems from the fact that as amr increases, the increasing contribution of kinetic energy 
over potential energy prevents gravitational collapse and thus inhibits star formation. As 
increasing a^r implies decreasing the clump's mass/size, as a^r oc L~'^ oc Mc (at fixed 
Mach number and density), star formation is thus highest in the largest and most massive 
clumps, which correspond to the smallest values of a^ir, in agreement with the conclusion 
reached in the previous subsection. Giant molecular clouds in the Milky Way, with masses in 
the range ~ 10^-10^ Mq, are generally weakly gravitationally bound structures, with virial 
parameters in the range a^r ~ 0.3-3, with slightly decreasing ayir with increasing cloud 
mass, as just noted (Heyer et al. 2009, Murray 2011). Accordingly, the rather low observed 
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SFR values might partly stem from the fact that most molecular clouds have in general 
virial parameter values slightly above unity, i.e. depart from perfect virialization, being held 
together partly by the confining ram pressure of turbulent flows in the ISM. 

(c) Thermodynamics of the gas 

For the non-isothermal case, we need to specify a temperature-density distribution (see 
eq. (18) of paper II). We choose no = 200 cm~'^ and Tq = 20 K at lO'' cm"^. Assuming that 
the velocity dispersion remains unchanged, we recalculate accordingly the speed of sound 
and the Mach number in the non-isothermal case (see eq. (20) of paper II). 

As seen in Fig. HJ taking into account the thermodynamics of the gas yields values 
of SFR'lr, for the same clump conditions, smaller by a factor of about ~ 1.5-2 than the 
isothermal case in the nearly constant SFR^rr region, which corresponds to small values of 
Omr, but leads to a significantly steeper decrease for increasing values of the virial parameter, 
in particular for Omr ^1- As explored in details in Paper II, this is a direct consequence 
of the density- dependence of the Mach number for non-isothermal gas. Indeed, when taking 
into account the thermodynamics of the gas, the clump is warmer at low densities, which 
implies lower Mach numbers than for the isothermal case (see Paper II), decreasing the SFR, 
as seen in the figure. 

As seen in the figure, for low Mach numbers, which correspond to small-size clouds, 
the value of i/cut has a significant impact on SFR'jj-, as already discussed in subsection 
(a), leading to substantial uncertainty on the SFR for such low- mass/small-size clumps. 
Fortunately, such clumps are predicted to contribute almost negligibly to star formation, as 
mentioned earlier. 

(d) The role of magnetic field 

The two bottom panels of Fig. H] show the SFR in the two magnetic cases corresponding to 
V^" = 1, 76 = 0.3 and V^ = 3, -yb = 0.1, respectively, adopting i/cut = 0.1. 

In order to decipher the contributions of the various physical effects, we have calculated 
the SFR in the magnetized case first keeping the same PDF as in the hydrodynamical case 
(as described by eq. H]) but then we have also taken into account the fact that the magnetic 
field can modify the width of the density PDF, as investigated for example in Molina et 
al. (2012). In a magnetized flow, the dependence upon the Alfvenic Mach number must 
be taken into account. Based on their numerical simulations, Molina et al. (2012) propose 
an analytical relation to predict the dependence of the variance of the lognormal density 
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Fig. 4. — Values of SFR^r^ as a function of a^r corresponding to various cloud parameters 
and A^ = 16 (top), 9 (middle) and 4 (bottom), in the isothermal, non-isothermal, hydro- 
dynamical and magnetized cases, for 2 values of ycut, as indicated on the figure, and our 
fiducial value b = 0.5. For the magnetized cases (right and left bottom panels), the solid 
lines correspond to the case where the magnetic field is not supposed to affect the density 
PDF, while the dotted lines display the results when the effect of the field on the width of 
the PDF is taken into account (see text). 
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distribution upon the rms Mach number in magnetized supersonic turbulent gas. The result 
depends on the index %. For % = 0, i.e. B independent of the density, they recover the 
hydrodynamical dependence (eq. H]). When 75 =0.5, i.e. B oc ^p, they find that, in the 
super-Alfvenic case, the variance is given by 



al = ln{l + b'M'^), (49) 

where P = Pth/Pmag = 2C2/K'- 

The basic effect of the magnetic field on the SFR, when assuming no variation of the 
density PDF compared with the hydro case, is illustrated by the solid lines in the bottom left 
and right panels, to be confronted with the SFR values displayed on the right third panel. 
As expected, the magnetic field tends to reduce the SFR by a factor of a few. The amplitude 
of this effect depends obviously on the field intensity, as explained in §2.3.31 For V^ = 1, 
7fe = 0.1 (not displayed here), the SFR is identical to the hydrodynamical case, while for 
V^ = 3, 76 = 0.1 (left bottom panel), it is reduced by a factor of less than 2 for A^ = 16 
and about 2 for A^ = 9. Stronger values of the magnetic field obviously lead to even lower 
SFR's. However, the values we used are within the range of the values typically measured 
in molecular clouds. 

The dotted lines of bottom left and right panels show the SFR when the infiuence of 
the magnetic field on the density PDF is taken into account, as described above. With our 
choice of parameters we find that at the mean density (3 = 0.036 for V^ = 3, % = 0.1 and 
P = 0.6 for V^ = 1, 7b = 0.3. The first value in particular significantly reduces the width 
of the density PDF. Clearly the impact of the magnetic field is much more pronounced and 
the SFR drops by a factor of several. We stress, however, that, as discussed in Mohna et 
al. (2012), the relation given by eq. fH^ only holds for 7^ = 0.5 and is is not working well 
at small and intermediate Mach numbers. Therefore, there is still considerable uncertainties 
here due to the lack of exact knowledge of the density PDF in the presence of magnetic field 
and also of the precise gas density-magnetic intensity relation. 

Finally, we also note that in the presence of a magnetic field, the coefficient (pt could 
possibly be further increased, reducing further the SFR. Indeed, it is likely that the magnetic 
field and the velocity field of the perturbation that eventually gives rise to the gravitation- 
ally unstable perturbation must be sufficiently aligned. Otherwise, after a weak density 
enhancement, the magnetic pressure will probably stop the contraction. 
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(e) The role of turbulence 

As in the theories of Kruinholz & McKee (2005), Padoan & Nordlund (2011) and Hopkins 
(2011), turbulence is at the heart of our star formation theory. We have shown in Pa- 
per I (§6.3) that, in the time-independent version of our IMF theory, turbulence, globally, 
has a negative impact on star formation. The stronger the turbulence, both at the global 
clump scale (as formalized by M.) and local (Jeans length) scale (as formalized by A^*), the 
smaller the global clump mass fraction encapsulated into gravitationally bound prestellar 
mass reservoirs. Indeed, although increasing the Mach number increases the width of the 
PDF, promoting the formation of small-scale structures, it leads to a decrease of the number 
of collapsing structures around the peak of the CMF (see papers I and II). As shown in the 
previous sections, however, star formation will keep processing during the clump's lifetime 
and time dependence changes the final clump's mass fraction ending up forming collapsing 
cores. The impact of turbulence on SFR^jrj is illustrated in FigEl This figure portrays 
the behaviour of SFR^rr at fixed cloud density (n=constant along each line) for different 
cloud sizes, as a function of Vq, thus of A^ (see eqn. flT6l) ). spanning a range A^ ^ 4-30. 
The time-independent results are displayed in Fig. [61 In this latter case, the SFR steadily 
decreases with increasing levels of turbulence, as mentioned above. As shown in Fig. 
however, including time- dependence in our formalism significantly changes this behaviour, 
with the SFR now increasing, although modestly, with the level of turbulence. The reason 
is twofold. First, increasing the Mach number shifts the peak of the CMF towards smaller 
masses (see Paper I), a behaviour exacerbated when time-dependence is considered, as dis- 
cussed in §2.2.41 Since small-scale structures have shorter free-fall times {tr oc p^^/^), this 
increases the number of small cores and thus globally increases the SFR (see eqn. fHO]) ). 
Second, massive mass reservoirs, considered to be stable in the static case, can now fragment 
into small-scale structures as their internal turbulent motions lead to new small-scale over- 
dense regions, globally increasing the number of collapsing cores. This positive impact of 
turbulence upon the SFR is in agreement with the results of PN but contrasts with the ones 
of KMK. As demonstrated in HCll, this decreasing dependence of the SFR with increasing 
Mach number in KMK theory (see their eq.(30)) stems from the missing pj^^^ term in their 
eq. (20), which stems from the fact that KMK assume that the critical free-fall time is equal 
to the free-fall time of the cloud. The dependence of the SFR upon b illustrated in Figl5]has 
already been addressed previously. 

The fact that, in the time-dependent formalism, turbulence is found to enhance star 
formation seems, at first glance, to be in contradiction with numerical simulations (see e.g. 
MacLow & Klessen 2004). However, as discussed in HC2011, the behaviour we infer for 
the SFR is qualitatively similar to the one inferred by Padoan & Nordlund (2011) (see also 
Federrath & Klessen 2012). In their simulations, these authors find that, for a given a^r, 
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the SFR increases with the Mach number. The critical point here is that when turbulence 
is too high (i.e. a^r is typically larger than 1), the kinetic energy will lead to an expansion 
of the clump, an effect which is neither included in the present analytical calculations nor 
in the turbulence in the box numerical simulations of Padoan & Nordlund (2011). We thus 
intuitively expect a non-monotonic behaviour of the SFR with turbulence. When ayir is 
small, turbulence enhances star formation because of the presence of converging motions that 
lead to density enhancements. On the other hand, if turbulence becomes too large {a^ir > 1 
or so) , it triggers a fast expansion of the cloud which reduces the density and quenches star 
formation, (note that Vazquez-Semadeni et al. 2003 do see such a non-monotonic behaviour 
with turbulence). 

(f) Summary of the results 

To conclude this section, we summarize the main results of our calculations at this stage. As 
mentioned above and seen in Fig. |3]and[5l for all typical clump conditions we have explored, 
which correspond to clump sizes ~0.5-20 pc and masses ~ 200-10^ Mq, A^ = 4-30, SFR^^j 
lies in the range ~ 0.01-5 in the non-isothermal case, even though the value of Hcut signifi- 
cantly affects these values at low Mach. Lower values of SFR^rr would imply either ycut of 
the order of 1/30 or even 1/100, which does not look very realistic, or very small clumps 
(<0.1 pc), for which the mass spectrum is drastically truncated and is completely incom- 
patible with the observed CMF/IMF (see Fig. 8 of paper II). According to our calculations, 
the value SFR^rr ^ 0.01 thus seems to be a lower limit for realistic clump conditions. As 
shown above, the larger SFR values are reached for the densest clumps and the largest/most 
massive clumps, which correspond to the smaller values of the virial parameter and are thus 
likely to be gravitationally bound structures. Significant star formation, however, is still 
predicted to occur in low-density clumps, providing they are large enough (i) to generate 
enough turbulence levels (high enough Mach numbers) and (ii) to largely exceed their typ- 
ical Jeans mass. We must recall, however, that SFR^rr is the SFR obtained assuming that 
self-gravitating fluctuations are replenished within one single free-fall time (0t = 1) and with 
a 100% efficiency of initial mass reservoir to star conversion during the gravitational collapse 
(e = 1). SFR^^f must thus be multiphed by a factor e/cpt < 1 to yield the effective value 
SFRff. The dependence of these two parameters upon the local physical conditions such 
as magnetic field intensity or radiative feedback remains very uncertain, but both processes 
tend to decrease the final star formation rate SFRfj, either by decreasing e or by increasing 
(pf As mentioned earlier, observational and numerical results suggest that e/(f)t ~0. 1-0.2. 
Therefore, according to our calculations, the star formation rate per free-fall time in dense 
molecular clumps with typical aforementioned properties is estimated to lie in the range 
SFRff ;^0. 001-1.0, for ycut = 0.1-0.3, for virial parameters in the range 0.3<a«r.<3, in 
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agreeinent with typical observed values (see e.g. Murray 2011 and §1] below). The cal- 
culations of a mean value of the SFR, integrating over the clump mass spectrum, will be 
addressed in 3S1 



4. Comparison with observations 

4.1. Star formation rate versus gas surface density 

As seen in the previous section, the SFR per free-fall time can vary by several orders of 
magnitude depending on the clump's conditions (mass, density, typical Mach number) but 
also, unsatisfactorily enough, on the values of ycut, b, (pt and e. These large uncertainties 
prevent a precise theoretical determination of the SFR in a cloud. Nevertheless, it is instruc- 
tive to compare our determinations with observational ones. Such a comparison has partly 
been done in HCll and will be extended here. The observed star formation rate per free-fall 
time in star forming molecular clouds in the Milky Way lies in the range 0.03 < SFRff < 0.3 
for clouds in the mass range 10^ <Mc/Mq< 10^, with Ai ~ 10-15 and a^r ~ 0.3-2, with a 
mean value (SFRff) ^ 0.16 (Murray 2011, Heyer et al. 2009). Looking at Galactic nearby 
molecular clouds and massive star-forming dense clumps, Heiderman et al. (2010) find SFR's 
in the range ~ 0.02-0.12 for molecular clouds and ~ 0.03-0.5 for dense clumps, yielding a 
mean value f« 0.1. These values are quite consistent with the ones derived above from our 
theory, while the typical rates obtained by Krumholz & McKee (2005) for similar conditions 
are about an order of magnitude smaller, as noted by Murray (2011) and Heiderman et al. 
(2010), a point we addressed in the previous section. 

Traditionally, the observationally determined star formation rate for a given clump (or 
region of a clump) of mass Mc, radius Re, lifetime Tq and mean dynamical time r9r is the 
projected SFR, i.e. the SFR per unit area : 

t, = ^ = SFE X (^) = SFRff X (^), (50) 

To T-Q Tff 

where Eg = Mc/nRl = {4:/3)pRc for a spherical clump is the clump gas surface density. The 
relationship between this projected SFR and the gas surface density is generally found to 
obey, within a large scatter, a power law relationship: 
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Extragalactic studies yield the well known Kennicutt-Schmidt relation, with Sioo = 0.16±0.6 
and A^ = 1.4 ±0.15. 
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Fig. 7. — Comparison of the SFR per unit area, E+, as a function of gas surface density, 
Eg, with various observational determinations: Heiderman et al. (2010) massive clumps 
(triangles) and molecular cloud YSOs (diamonds+squares), Gutermuth et al. (2011) (areas 
bracketed by dashed lines). The four solid lines correspond to four clump sizes, namely 
i?c = 0.5, 2, 5 and 20 pc (right to left). 
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Figure [7] displays S^ predicted by eqn. f H2|) . as a function of surface density, S^, for 
four typical clump sizes, namely Re =0.5, 2, 5 and 20 pc, both in the isothermal and non- 
isothermal case. The impact of the uncertainties due to i/cut is shown by the two calculations 
with i/cut = 0.1 and 0.25, respectively. The clumps are assumed to follow Larson's (1981) 
relations and thus have velocities and densities given by eqn. (135|) and (!36l) . with mq = 1 and 
do = 0.1-5, i.e. n ~ 10-10^5- The corresponding global and local Mach numbers entering our 
theory, Ai and A^*, are consistently derived from these values. We have taken e/0( = 0.1 
in the figure. Also displayed on the figure are the data of Heiderman et al. (2010), both 
for the class I and flat-SED YSOs and the massive clumps as well as the location of the 
observations of Gutermuth et al. (2011), as bracketed by the two large diamond areas. 
These latter authors determined the SFR in eight nearby molecular clouds, sampling gas 
surface density regions, inferred from near-IR extinction mapping, ranging from ~ 15 to 
300 MqPc"^. Their analysis is similar to the one by Heiderman et al. (2010) but, as they 
examined larger clouds, put stronger statistical constraints on the SFR vs gas surface density 
determinations than that study. Moreover, these authors include pre-main-sequence Class 
II objects in addition to Class 0, I protostars, from Spitzer and 2MASS surveys. 

Several conclusions can be drawn from the theory-observation comparison displayed in Fig. 

m 

(i) all calculations reproduce well the observed values, with the scatter, and exhibit a 
direct correlation between the density of star formation and the gas density, in agreement 
with the observational determinations. 

Clearly, it seems difficult to provide a universal relation such as eq. (l5T|) over the whole 
density range, given the strong dependence, both for the slope and the normalization, of the 
SFR upon the properties of the clump, i.e. its mass and the nature of turbulence forcing. 
What seems to be more robust is that, at very high density, S^ > 300 Mq pc~^, the theoretical 
relations seem to merge towards a strongly super-linear behavior with an exponent A^ ~ 2 
(an exact value is difficult to infer as it depends on the clump size) and a normalization 
S300 ~ 32. Star formation is thus predicted to increase basically quadratically with density 
and thus to be largely dominated by the contribution of the densest (likely to be bound) 
clumps. Given the steeply decreasing number of clumps with increasing density, however, 
such dense clumps represent only a modest ( < 20%) mass fraction of the clouds (Kainulainen 
et al. 2011). This slope becomes progressively steeper as density decreases, the steeper the 
smaller the clump, as expected from the analysis conducted in the previous section. 

This A^ ~ 2 value for the slope is in excellent agreement with the values A^ = 1.67 to 
2.67 determined by Gutermuth et al. (2011) in eight nearby molecular clouds. Lower values, 
including Kennicutt-Schmidt like values, A^ ~ 1.4, with a a star formation threshold are 
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excluded by these authors. 

(ii) the dispersion in the observed SFR determinations is well explained by variations 
among the clump sizes, not mentioning possible variations of the (e/0t) efficiency factor and 
of the contribution of different turbulence modes, as formalized by the factor b, between 
clumps. 

(iii) the theory, as the observations, predict a severely decreasing SFR below a typical 
value Sjj ^ 100 Mq pc~^, i.e. Nh^ ~ 5 x 10^^ cm"^, n ^ 2500 {R/l pc)"^§, which corresponds 
to a visual extinction Ay ~ 7 (Draine 2003), except for the largest (20 pc) clumps. We recall 
that such a drop in the SFR naturally arises from our theory, without invoking any density 
threshold. 

(iv) there is no real density threshold in star formation. Indeed, according to our calcu- 
lations, star formation is predicted to keep processing, although with a steadily decreasing 
rate, even at densities lower than the aforementioned value, but only, at non negligible rates, 
in the largest clumps. This was expected from the discussion in §3. 2. 2( a). In this low-density 
domain, our values can be compared with the observations of Gutermuth et al. (2011). As 
mentioned earlier, these authors examined larger clumps than previous studies and found 
larger SFR values than these latter at low density. The eight molecular cloud sample ex- 
amined by these authors covers a range of clump sizes and surface densities Lc ~ 4-20 pc, 
Tjg ~ 22-71 Mq pc^^. In this range, they determine S^/S^ in the range (3 x 10~'^)-(5 x 10^^) 
pc^ Mq^. These values can be directly confronted with the ones inferred from the two upper 
solid lines, which correspond respectively to -Re = 5 and 20 pc, in the appropriate surface 
density range. This yields S^/S^ ^ 2.0 x (10~^-10~^) for the non-isothermal case and val- 
ues larger by about a factor 10 for the isothermal case, for our fiducial value of turbulence 
forcing, b = 0.5. Remembering that both observational and theoretical determinations are 
subject to significant uncertainties, the agreement can be considered as satisfactory. 

This reinforces our suggestion that the star formation rate smoothly decreases with 
decreasing gas density and that there is no real density threshold for star formation. The 
severe drop in the SFR, which is predicted to occur around Eg ^ 100 Mq pc~^ for ~pc-size 
clumps, simply reflects the basic mechanisms mentioned in §3.2.21 star formation can still 
occur, although at significantly lower rates, in low-density clumps providing these latter, as- 
suming they follow Larson's relation, are large enough to (i) significantly exceed their typical 
Jeans length and (ii) generate enough large-scale turbulence, i.e. have large enough Mach 
numbers, to generate a large enough spectrum of density fluctuations, and thus of prestellar 
cores. According to the observed clump size distribution, dN^/dLc oc L^^'^ (Kainulainen et 
al. 2011), however, such clumps are very rare. This stresses the need to observe very large 
areas at low-density in order to get statistically significant SFR detections. 
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A threshold for core formation was also ruled out in the SCUBA survey of Perseus by 
Hatchell et al. (2005). These authors stress the existence of submm cores, identified from 
IRAS as likely Class I objects, down to Ay = 3, i.e. Nh2 ~ 10^^ cm~^, and demonstrate 
the steeply decreasing probability of finding a core with decreasing column density. This is 
also consistent with the analysis of Andre (2012, Fig. 3), who find a drastically decreasing 
but non-zero number of class-0 prestellar core candidates below about the aforementioned 
Ay ~ 7 value. 

In contrast, several studies of various local clouds (Onishi et al. 1998, Johnstone et al. 
2004, Enoch et al. 2007, Lada et al. 2010, Andre et al. 2010) do find a steep break in the SFR- 
gas density relationship around the aforementioned ~ lOOMopc^^ value, that they identify 
as a density threshold for star formation. It should be stressed, however, that all these studies 
still find pre-stellar cores below this limit, suggesting that star formation continues, although 
with a much smaller efficiency, down to the limit of the surveys, i.e. "Eg ^ 10 MqPc^^. As 
mentioned above, a noticeable exception is the recent work by Gutermuth et al. (2011), 
who explore larger areas and who include pre-main-sequence Class II objects in addition to 
Class 0, I protostars, from Spitzer and 2MASS surveys, and who do not find evidence for a 
column density threshold for star formation from 15 to 300 Mq pc^^. These authors explain 
the different result by the very small number of Class 0, I protostars at low gas column 
densities, and thus the low statistics in previous surveys, including the Spitzer-derived data 
of Heiderman et al. (2010). Interestingly, by comparing the SFR derived from YSO surface 
densities and the ones derived from far-IR luminosity, as commonly used as a proxy of the 
SFR, they find that the LpiR-derived SFRs are systematically about an order of magnitude 
below the YSO-derived ones. This is confirmed for instance by the recent analysis of W40 
and Serpens South with Herschel by Maury et al. (2011), who have access to Class 
objects and find SFRs an order of magnitude higher than the typical SFRs observed for 
embedded infrared clusters. This suggests that, if anything, SFR determinations based only 
on Class or Class I objects underestimate the real SFR by a fair amount. Therefore, the 
aforementioned identification of a cut-off in the SFR below ~ 100 M© pc~^ might in fact 
be due to the very small number of Class 0, I protostars at low gas column densities and 
the too limited statistics (too small fields of view) in the surveys. This is supported by the 
recent analysis of Bressert et al. (2010) who show that using only near-infrared detections 
and small fields of views allows identifications of young stars only in the densest parts of the 
clouds, a bias which can be corrected with Spitzer data and larger surveys. 
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4.2. Molecular cloud vs galactic SFR 

As already found in the observations of Heiderman et al. (2010) and Evans et al. (2009), 
our SFR-gas relation, characteristic of Galaxy star-forming molecular clouds, lies above the 
Kennicutt-Schmidt relation by about a factor of 10 at high density. It should be remembered, 
however, that the Kennicutt-Schmidt relation derives from extragalactic determinations and 
is averaged over much larger regions (~ kpc-size) than individual molecular cloud complexes 
( < 10-100 pc). Such an average determination includes star-forming regions but also diffuse 
molecular gas or atomic gas that is not forming stars, leading to overestimate the amount 
of gas counted as star-forming gas (Heiderman et al. 2010). Recent SFR determinations in 
nearby galaxies indeed show that the measurement size scale used changes significantly the 
SFR-gas surface densities relation (Liu et al. 2011). Moreover, SFR in CMC's are determined 
from inventoring the YSOs and assuming a star formation timescale of ~Myr while galaxy- 
averaged SFR's are derived from conversion of a FIR flux into a mass growth rate assuming 
a timescale of 10 to 100 Myr, naturally leading to a smaller SFR for a given average gas 
density. Since, as shown in the present study, the SFR depends on the clump's mean density 
and significantly decreases for low-density clumps, this provides a plausible explanation for 
the difference between CMCs and extragalactic SFR-gas quantitative determinations. In 
fact, a correct determination of the SFR-gas relation at extragalactic scales should take into 
account the density-dependence of the SFR-gas relation. The same explanation applies to 
the characteristic value of the SFR inferred for the Milky Way, found to be about an order 
of magnitude smaller than in CMCs. A likely explanation is that most of the mass of the 
clouds which compose the Calaxy is at low column density, with only a very small fraction in 
condense structures at high enough density for efficient star formation (see e.g. Kainulainen 
et al. 2011). Indeed, large-scale clouds which compose the Calaxy have densities around 
~ 10^^ cm~^ (n < 100 cm~'^). Therefore, most of the clouds in the Milky Way no longer 
produce stars. On the other hand, our calculations suggest that the power dependence of 
the SFR upon the gas surface density in the (high-density) efficient star formation regime is 
not drastically different for Calactic and extragalactic determinations, since our exponent at 
high density is A^ ~ 2 whereas the one characteristic of the Kennicutt-Schmidt relation is. 
A^ = 1.4±0.15 (see Fig. [7]). Again, such a universal SFR-gas dependence naturally emerges 
from a picture of star formation being initiated by turbulence-induced fluctuations and the 
star formation rate being determined by the free-fall time of these fluctuations - not of the 
clump ! - once gravity takes over, leading to a strong dependence of the SFR upon the gas 
density. 
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4.3. Other dependences of the star formation rate 

In the previous section, only the SFR dependence upon the column density was exam- 
ined. In the present section, we consider the SFR dependence upon the clump's mass and 
upon the column and volume densities over free-fall time. The results are compared with 
the sample of Heiderman et al. (2010). 



4-3.1. Star formation rate versus clump mass 

Figure [Uportrays the star formation rate, M^, as a function of the cloud mass, Mc, for the 
isothermal and non-isothermal case. Each line corresponds to the same clump sizes as in the 
previous figures, i.e. Re = 0.5, 2, 5 and 20 pc, with the normalization density do increasing 
along the line from do = 0.1 to 4.5. We note the large scatter in the predicted M^ for a given 
clump mass, in particular for the non-isothermal case. The data of Heiderman et al. (2010) 
are displayed in the figure for comparison. The theoretical determinations agree very well 
with the data, indicating that, within our theory, the range of clump characteristics we have 
investigated reproduce well the observationally determined SFR values with the observed 
scatter. The dashed line in th figure displays the relation derived by Lada et al. (2012), 
M^ = (4.6 X 10~^) (^^) M0yr~^. The present calculations show that this relation, although 
broadly correct, may not be an accurate representation of the real star formation process, 
as this latter not only depends much more steeply on the clump's mass but also strongly 
depends on the clump's size. 



4.3.2. Dependence of star formation rate on free-fall time 

The values of star formation rates per mean free fall time for a sample of clumps can 
be directly inferred by determining the slope of the relation S^ vs (Sg/r^j) (eq. flSQj) ) or, 
equivalently, of the volume density relation n^ vs (n^/r^j). A universal value of the SFR 
would correspond to one single, constant value of the slope for the whole sample. 

To confront such observational determinations with our theory, we first vary the clump's 
size for a fixed density normalization parameter do. Indeed, assuming Larson's relation for 
our clumps, i.e. p oc Rc^'^, yields S oc R^'^ and r// oc R^'^^ for the column density and 
free-fall time, respectively. Therefore the column density over free-fall time, ^/Tff is almost 
independent of the clump's size. The corresponding S^ vs (Tig/r^j) curves are thus almost 
parallel to the y-axis, making the dependence upon do and Re straightforward. Seven values 
of do are computed, namely 0.5, 1, 2, 3, 4, 5, 10 as well as eleven sizes Re =0.5, 1, 2, 4, 6, 
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Fig. 8. — Star formation rate, M^, as a function of the cloud mass, Mc, in the isothermal and 
non-isothermal case. Each line corresponds to clump sizes Re = 0.5, 2, 5 and 20 pc, from 
left to right, with the normalization density d^ increasing along the line from do = 0.1 to 4.5. 
The squares, triangles and diamonds represent the data from Heiderman et al. (2010). The 
dashed line displays the fit of Lada et al. (2012). 



-42- 






non — isothermal, y„„t=0.25, b = 

3fe \' 



■0.5 




12 3 4 

log(I/U) (M. yr- kpc-^) 

non — isothermal, y<,ut=0.1, b = 0.5 



ex 




m 
o 


-2 
C 


1 

) 1 2 




3 






log(I/t 


„) (M. 


yr" 


kpc 



Fig. 9. — SFR per unit area, S^, as a function of (Sg/r°j) for 7 values of do =0.5, 1, 2, 3, 4, 
5, 10 (left to right) and 11 clump sizes Re =0.5, 1, 2, 4, 6, 8, 10, 12, 14, 16 and 18 pc (bottom 
to top on each line, as indicated by the crosses). Each curve corresponds to a constant do. 
The highest values of do and Re correspond to the highest SFR's. 
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8, 10, 12, 14, 16 and 18 pc. Figure |9] displays the theoretical predictions for S^ vs (Sg/r^j) 
as well as the observational determinations of Heiderman et al. (2010). As obvious in the 
figure, and as expected from the previous subsection, the overall agreement between our 
predictions and the data is very good and the observational scatter is well reproduced by 
the explored variations of clump sizes and densities. Note that, given the large spread in the 
data and the uncertainties discussed previously in the theory, it is difficult to assess more 
than a qualitative agreement. 
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Fig. 10. — Sj, as a function of the clump radius, r. There is a clear trend of S^ strongly 
increasing with decreasing clump's size, particularly if the massive clump determinations 
(triangles), which tend to differ from molecular clouds YSOs (diamonds and squares), are 
excluded. 

Interestingly, the confrontation between our results and the observed SFR values from 
Fig. |9] leads to the conclusion that the size of star forming clumps should decrease with 
increasing column density over free-fall time. Typically, the clump's size should be less than 
about 2 pc for log(S/rjj) ~ 3 and less than about 0.5 pc for log(S/rjj) ~ 3.5 (with a 
dispersion of a factor 2-3). To verify whether this trend is indeed real, we have plotted S^, 
as a function of the clump's radius (Fig. [10]). Indeed, Fig. [10] shows a clear trend for the 
highest values of S^, been obtained for the smallest clouds. The reason of this behaviour is not 
completely clear but appears to be an important element to understand the SFR dependence. 
A likely explanation is that large clouds never become dense enough, i.e. massive enough 
(see e.g. Fig. 8). More precisely, as massive clouds collapse they undergo local gravitational 
fragmentation in dense clumps where most of the star formation occurs. 
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Fig. 11. — Volumic SFR density, rij,, as a function of {ug/r^j-) 
Figini Each line represents a fixed value of c/q. 



for the same conditions as in 



Finally, for completeness we also display n^ vs {ng/r^p) in Figure [TT] for the same cloud 



parameters as in Fig. IH] (the various curves correspond to a fixed do) while Fig. [T^ shows the 
results for the same cloud parameters as in FigJTl For the data of Heiderman et al. (2010), 
we derived the volume density from the surveyed area as (assuming spherical clumps of area 

According to the present analysis, there is clearly no "universal" slope value neither 
for the S^ vs (Sg/r9r) nor for n^ vs (ng/r^f) relation. Instead, these relations significantly 
vary with the clump's size and density. Indeed, below ^/Tfj < 1000 M0yr~^kpc~^ there is 
between 1 and 2 orders of magnitude variations in the SFR. Even in dense regions, corre- 
sponding to n/Tff > 100 cm~'^Myr~^, the mean value of the SFR, as measured by the slope 
of the illustrated relations, typically varies between ~ 3% and > 10%. This analysis is in 
contrast with the claim by Krumholz et al. (2012) for a universal SFR value, of about 1% 
from clouds to starburst galaxies. As mentioned earlier, molecular clouds are very clumpy 
structures with orders of magnitude of density variations. Therefore, it seems difficult to 
invoke an average timescale, defined at the cloud's mean density, as the relevant timescale 
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for star formation. 
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Fig. 12. — Volumic SFR density, n^, as a function of {ug/r^j-) for the same conditions as in 
FigJTl The four sohd hues correspond to four clump sizes, namely Re = 0.5, 2, 5 and 20 pc 
(right to left). 



4.4. Discussion: The case of filaments 



According to various recent studies, in particular the ones conducted with Herschel, 
star formation is found to occur in filaments (Andre et al. 2010, Men'shchikov et al. 2010, 
Ward-Thompson et al. 2010, Molinari et al. 2010, Hill et al. 2011). It is well known that 
(unbound or bound) molecular clouds exhibit a filamentary structure, as expected from both 
hydrodynamical and MHD compressible (shock dominated) turbulence even in the absence of 
gravity (Padoan et al. 2001, Nakamura and Li 2008, Federrath et al. 2010, see also Banerjee 
et al. 2009). It is thus not surprising that the locations of star forming populations tend 
to follow this morphology. We recall that in our formalism, the prestellar cores, progenitors 
of the protostars identified as Class objects, take birth in turbulence induced clumps of 
enhanced density, and isolate themselves from the surrounding medium under the action of 
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gravity (Hennebelle & Chabrier 2008, Chabrier & Hennebelle 2011). For sake of simplicity, 
the clumps and the cores in our theory are assumed to have a spherical geometry. Non-linear 
simulations of the collapse of filaments show that the resulting fragments are nearly spherical 
(Inutsuka and Miyama 1997). It thus seems reasonable to assume spherical collapse for the 
cores. The reality might be more complex for the star-forming clumps, which may have 
a flattened, filament-like geometry, some of these filamentary clumps becoming themselves 
gravitationally unstable as they both accrete mass and dissipate turbulent energy, yielding 
eventually global collapse of the clump. Interestingly, the very same Herschel observations 
reveal that protostellar cores seem to be present only in the gravitationally bound clumps 
(Andre et al. 2010, Molinari et al. 2010, Arzoumanian et al. 2011), suggesting that star 
formation takes place dominantly in such clumps, prone to gravitational fragmentation. 
Gravitational contraction will thus cause filamentary clumps to fragment along the axis. 

This possible peculiar role of filaments in star formation might affect our results by 
numerical factors of order of a few but it does not fundamentally modify our theory of 
clump and prestellar core formation, and of the resulting CMFlj . In our theory, the (scale- 
free) criterion for clump formation is simply set up by an arbitrary critical density threshold 
6c = log(pc/p) (see §3 of paper I). The argument above thus suggests to choose for this 
threshold condition for the star-forming clumps - identified as the aforementioned filaments 
- the density at which these latter become gravitationally unstable and thus exceed a critical 
mass per unit length Merit = 2(o")^/G', where (a) includes both the thermal and non-thermal 
velocity dispersion contributions (Ostriker 1964, Larson 1985, Inutsuka & Miyama 1992, 
Fiege & Pudritz 2000). For a temperature T ~ 10 K, this corresponds to a mean column 
density Nh2 — 10^^ cm~^, similar to the value identified in §4.1 for efficient star formation 
(point (iii)), about 5-10 times the typical mean density of the surrounding cloud. This 
particular choice for the clump formation in our theory will not affect the resulting core 
mass function (CMF), determined by the virial condition for spherical collapse, as mentioned 
above. 

Furthermore, a closer inspection of Table 1 and Fig. 8 of IIC09 (see also Fig. 2) shows 
that small-size, low-density (< 3000^) clumps, for which the CMF is very narrow and thus 
star formation is neghgible, are unbound or only marginally bound ((-Efcm + Efh) ^ \Egrav\), 
whereas the ones above about this density are bound and thus prone to collapse. Therefore, 
our theory suggests that the clumps in which star formation is taking place dominantly 
should be denser than n ~ a few 10'^ ^, i.e. Nh2 — 10^^ cm~^ for pc-size clumps, in fairly good 
agreement with the aforementioned critical value inferred for filaments. This is consistent 



^The characteristic Jeans mass for a nearly isothermal filament (7^1) differs only by a factor 0.6 from 
the one for a sphere (Larson 2003) 
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with the results found in §3.2.2[ that show that the SFR is highest in the largest and densest 
clumps, with a^r ^ 1- Note that, given the scale-free nature of the turbulent flow responsible 
for the clump formation, such a minimum density for efficient star formation must apply at 
all scales and thus must vary accordingly with the clump's size. The aforementioned value 
typically applies for a 1 pc size clump and should increase as ~ -R^°'^ for smaller sizes, 
providing Larson's relations apply. 

Therefore, a picture where star formation occurs dominantly in relatively dense, rather 
bound clumps, with mean density n > 3000 5, in agreement with observational results, nat- 
urally emerges from our theory. The ffiamentary nature of the clumps does not significantly 
affect the results, except for possible numerical factors of order unity. 



5. Influence of the clump distribution and evolution 

So far, we have considered the core mass spectrum produced within a single clump with 
fixed physical quantities (mass, size, density). In real situations, there will be a distribution 
of evolving clumps. In this section, we investigate the impact on the gravitationally bound 
core mass spectrum of the time dependence of the cloud parameters as well as having a 
clump distribution rather than a single clump. In order to distinguish the influence of these 
two effects, we first consider the time evolution of a single clump before investigating the 
effect of a distribution. 



5.1. A simple model of clump evolution 

One of the intrinsic difficulties in inferring the mass spectrum of the self-gravitating 
fluctuations that may eventually lead to stars is the choice of clump parameters (density, size, 
velocity dispersion...). As seen from our calculations, these clump properties have a strong 
impact on the mass spectrum. The peak position, for instance, is determined by the clumps 
Jeans mass, thus mean density, temperature and Mach number. Clumps, however, are not 
well deflned entities, in particular because they constantly grow in mass by accretion, so that 
their density and size evolve with time. Under these conditions, it is clear that the choice 
of a flxed set of clump parameters is a simpliflcation. A proper investigation of the clump 
evolution would require to develop a fully self-consistent model for clump formation and 
evolution, including growth by accretion, gravitational contraction or turbulence dispersion 
as well as spatial inhomogeneity. There is no such model of clump evolution yet; we thus 
consider a very simple prescription, essentially inspired by observations. 
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Fig. 13. — Clump evolution for M = 10^^ Mq yr^^ and e/0i = 0.1. Top panel: total cloud 
mass (solid line), gas mass (dotted line), star mass (dashed line). Middle panel: Cloud radius 
over Jeans length. Bottom panel: mass spectrum after different times: t =1, 4 and 7 Myr 
(dotted lines from bottom to top) and t = 10 Myr (solid line). The dashed line displays the 
SCIMF 
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Fig. 14.— Same as Fig. [15] for M = 10"^ Mq yr"^ and e/0i = 0.3. 
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Fig. 15.— Same as Fig. [15] for M = 10"^ Mq yr"^ and e/0i = 0.1. 
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Fig. 16.— Same as Fig. US] for M = 10"^ Mq yr" 
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As mentioned above, instead of setting a density and size, we assume that the clump 
is accreting at a fixed accretion rate, M. Starting with a mass M{t = 0) = Mo, the clump 
mass grows as M^. = Mq + Mt. This mass entails the mass of the gas, Mg, and the mass 
of the stars, M^, that formed in the clump, M^ = M^ — Mg. We further assume that the 
clump follows Larson- type relations as stated by eqs. (135!) with do = 3 and uq = 1, which 
correspond to the canonical values of these relations. Then, knowing Mg{t), we can infer the 
gas mean density, ng{t) and thus the cloud radius Rdf). Since at a given time, t, the clump's 
parameters are well determined, we can calculate the SFR from our CMF/IMF theory, as 
stated by eqs. P2l) - fj43l) . We can thus estimate the gas mass fraction converted into stars 
between t and t + dt as: dM^:{t) = Mg x SFRff x dt. The total mass of stars at time t, 
M^{t), is thus 

M,(t) = f Mg{t)SFRff{t)dt. (52) 

Jo 

We typically integrate from to 10 Myr with time steps dt = 10^ yr. 

The accretion rate onto the molecular clumps is a difficult quantity to infer observation- 
ally and thus is not well known. Using observations of the LMC, Kawamura et al. (2009) 
and Fukui et al. (2009) propose that giant molecular clouds with a mass of ~a few 10^ Mq 
are accreting at a rate M =(l-5)xl0~^ Mq yr~^. To estimate the accretion rate on smaller 
clumps, we simply use the Larson relations, 

M M 

M ~ — ~ 

Tc RJia/Vs) 

, 1.3+17 Uq j Ait , \ 2-^ ,^ , 0.7-317 



M^^-^(— rfo) (Ipc)^^ (53) 



1.3+17 



o / M \ 2.3 cr X / n \ ^ _ 1 

"- '" [wm;) 'Oki^HiT^) *^«="- • 

where the numerical constant has been estimated for rj = 0.4. For a GMC of mass 10^ Mq, 
this estimate yields an accretion rate M ~ 5 x 10~^ Mq yr~^, slightly lower but nevertheless 
close to the estimate of Fukui et al. (2009). For a cloud of 10^ M©, we get 2 x 10^"^ M© yr"\ 
We have thus explored three typical values for the accretion rate, namely M = 10~^, 10^'^ 
and 10~^ Mq yr~^, which correspond to clumps of approximately 10^, 10^ and 10^ Mq. 

Figure fT3] portrays the results for M = 10^'^ Mq yr^^ and e/0t = 0.1. The top panel 
shows the total mass of the clump (solid line), the mass of gas (dotted line) and the stellar 
mass (dashed line). While the total mass grows linearly with time, the mass of gas, which 
was initially almost equal to the total mass, increases less and less rapidly as star formation 
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is occurring. At the beginning of the process, the mass of the cloud is small, and the amount 
of stars formed remains fairly limited. After t ~ 4 Myr the stellar mass starts increasing 
significantly and after t ~ 8 Myr, the masses of gas and stars are roughly equal. After this 
time, the gas mass basically saturates, as the gas added to the cloud is rapidly converted 
into stars. As a consequence, the size of the cloud (middle panel) increases rapidly at the 
beginning and much more slowly as star formation proceeds. The bottom panel displays 
the stellar mass spectrum at different times, namely after 1, 4 and 7 Myr (dotted lines from 
bottom to top) and after t = 10 Myr (solid line); the dashed line correspond to the SCIMF. 
As seen in the figure, the shape of the mass spectrum remains about the same over time but 
of course the integral, i.e. the total number of stars, strongly increases, a consequence of the 
growing size of the clump, as mentioned above. 

Figure HM displays similar results for the same accretion rate, M = 10~^ Mq yr^^, but 
e/0f = 0.3. The behaviour remains qualitatively similar except that, as expected, stars form 
faster, since the SFR is three times larger (see eqn. ( H3|) ). and thus the gas mass remains 
smaller by a factor of about 2 compared with the previous case. Consequently, the size of 
the clump is slightly smaller and less small-mass objects form because the Mach number 
{Ai oc K^) is slightly smaller (see Paper I). 

Figure [15] displays the results for M = 10~^ Mq yr~^. The SFR is slightly smaller 
than for M = 10"'^ M0 yr^^ because the clump grows more rapidly and thus is less dense. 
Interestingly, the mass spectrum remains very similar, except for the slightly larger number 
of low-mass objects, a consequence of the larger Mach number. Indeed, as discussed in paper 
I (eq.(47) of paper I), there is a partial compensation between the Jeans mass and the Mach 
number scale dependences as a function of the cloud's size/mass so that the location of the 
peak of the CMF remains almost unchanged. It is interesting to compare the top panel of 
this figure with the top panel of Fig. 5 of Vazquez- Semadeni et al. (2007), which shows the 
total gas+stellar masses in a simulation of colliding flows. The general behaviour and even 
some of the details of the simulation are very similar to the results portrayed in Fig. [151 (i) 
the total mass of the cloud in both cases is a few 10^ M©, (ii) the time at which stars and gas 
masses are equal agree within a factor of about 2 (particularly if we choose for the beginning 
of the cloud formation in the simulation the time t ~ 10 Myr, which seems more accurate 
than t = 0, (iii) when stars start forming efficiently, the mass of the gas remains nearly 
constant both in the simulation and in the present model. The main differences appear in 
the evolution of the stellar mass, which is more sudden in the simulations and faster at the 
beginning than in the later phases, in opposite to the present results. Moreover, while in the 
model the mass of the gas always increases, it decreases in the simulation. This is likely due 
to gravity whose effects on the cloud dynamics are ignored in our model. 
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The last case we investigate is presented in Fig. [16] which displays results for M = 10~^ 
Mq yr~^. The mass of the clump is only about 400 Mq. The mass spectrum peaks at about 
the same mass, for the reason mentioned above but the mass spectrum is much narrower 
with a significant deficit of both massive and low-mass objects. This stems from the fact 
that for such small clumps, the Jeans length becomes comparable to the size of the cloud, as 
discussed earlier, drastically decreasing star formation (see e.g. Fig. 2, and Fig. 8 of HC09). 



5.2. A distribution of time-dependent star forming clumps 

It is now relatively well established that not only CO clumps (Heithausen et al. 1998, 
Kramer et al. 1998) but also infrared dark clouds (Peretto & Fuller 2010) follow a mass 
distribution A/''^' = dN'^^ /dM oc M~'^''^ with 'jd ~ 1.7. It is important to recall that our 
theory predicts that this exponent, as well as the one for gravitationally bound cores when a 
gravitational collapse condition is taken into account, is a direct consequence of turbulence 
and is indeed related to the index n' of the powerspectrum of logp by the relation 'y^i = 
1 + a; = 3 — n' /3 (see paper I). Therefore, the present theory naturally takes into account 
the clumpy structure of the gas and in principle no further calculations are needed. In other 
words, for spatially close enough clumps surrounded by sufficiently dense gas, there is no 
need to sum up over a clump population. While strictly true as long as i) the gas is barotropic 
and ii) the density PDF is lognormal, this is no longer the case in real situations. Indeed, the 
diffuse ISM does not follow our eq. ( 13T|) . as the atomic hydrogen which fills up the Galaxy 
can be as warm as 10^ K. Therefore, for clumps which are spatially well separated, the in- 
between gas is not isothermal and we must sum up over a clump population. The scale at 
which this happens is not well determined. Numerical simulations by Vazquez- Semadeni et 
al. (2007), Hennebelle et al. (2008), Heitsch et al. (2008), Banerjee et al. (2009) suggest 
that the warm gas is tightly mixed with the cold gas at scales of about a few parsecs. 

We calculate the mass spectrum resulting from the distribution of the time- dependent 
clumps described in the previous section. For that purpose, we proceed as follows. We 
perform a series of time- dependent clump calculations, choosing values of M between 10~^ 
and 10~^ Mq yr~^ using logarithmic intervals. We stop the integration when M* = Mg in the 
clump and we adopt for the mass of the clump the value M^ = M^. + Mg. Then we sum up 
the corresponding mass spectra after multiplying by the aforementioned clump distribution 
M'^\Md) oc M'^^K The final core mass spectrum is given by 



J^tot= MMlMa)VdM''\Ma)dM,u (54) 

Jo 

where MmX-^ci) = dn/dMc is the mass spectrum of self-gravitating fiuctuations (i.e. cores). 
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Fig. 17. — Integrated mass spectrum of self-gravitating objects (prestellar cores) for three 
values of 'jci, the index of the mass spectrum of the star forming clumps. Top panel: 7c/ = 2, 
middle panel: 'jci = 1-7, bottom panel: '-fci = 1-5. The three solid curves correspond to three 
different values for the upper value of the accretion rate, namely (from the leftmost line to 
the rightmost one): Msup = 10"\ 10"^ and 10"^ M© yr^^ The dash-line is the SCIMF. 
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i.e. the CMF, for a clump of mass Md and volume Vd, as given by eqn.(2) (see HC08 and 
HC09) and Mg^p is the largest mass in our sample which corresponds to the largest accretion 
rate. 

At this stage, it is worth stressing that, since 'jd < 2, most of the mass is contained 
in the most massive clumps; 'jd = 2 then corresponds to a critical case. We have thus 
considered three values of jd, namely 1.5, 1.7 and 2. Figure [T71 displays the results. 

For 7ci = 1.5 and 'jd = 1-7, the mass spectra are pretty similar to the results obtained 
in the previous section for a single clump. The reason is that most of the mass is contained 
in the most massive clumps. Interestingly, varying the upper value of the accretion rate by 
two orders of magnitude has only a modest impact on the upper mass part of the total mass 
spectrum while below about the mean Jeans mass, the low-mass part of the distribution 
quantitatively varies by orders of magnitude, the number of low-mass objects, in particular 
brown dwarfs strongly decreasing with decreasing accretion rate. Values of Msup < 10~^ Mq 
yr~^ yield a strong deficit of brown dwarfs compared with the SCIMF, representative of the 
observed population. As mentioned above, this stems from the smaller cloud size, thus the 
smaller Mach number with decreasing accretion rate, yielding a lack of overdense small-scale 
structures, progenitors of the brown dwarfs (see HC08 and HC09). 

For 7c; = 2, the difference with the mass spectrum obtained for a single clump is more 
pronounced. In particular at intermediate mass, the distribution tends to deviate from the 
SCIMF. This is due to the mass of the gas being equally distributed between small and 
large clumps. As the CMF which results from small clumps peaks at larger masses than the 
CMF produced in massive clumps, this tends to create a small bump at intermediate masses 
(M ~ IOMq). 



6. Conclusion. A paradigm for star formation and the initial mass function 

6.1. Summary of the results 

In this paper, we have generalized our previous formalism and developed an analyti- 
cal time- dependent theory of the stellar initial mass function. Although not very different 
from the time-independent prestellar mass spectra for given sets of cloud initial conditions 
(temperature, density, Mach number), the time-dependent ones are not quite identical to 
the former ones. The most significant differences are a slightly steeper slope at large masses, 
which arises from the time dependence of the characteristic Jeans mass, and the position 
of the peak of the CMF which is shifted towards smaller masses by a factor of the order of 
2 to 3, depending on the initial conditions. This is a direct consequence of the small-scale 
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structures being rejuvenated several times during the collapsing process of the larger ones. 
For all explored typical clump conditions, the resulting CMF, while slightly too narrow when 
considering an isothermal equation of state, agrees very well with the shifted Chabrier system 
IMF (i.e. the Chabrier IMF shifted by a factor of ~ 3 to take into account the core-to-star 
mass conversion efficiency during the collapse) when including the thermodynamics of the 
gas, confirming the results obtained for the time-independent theory. 

This time-dependent theory enables us to derive an expression for the star formation rate 
for molecular cloud conditions. Confrontations with previously published SFR theories show 
some similar trends but quantitative or even qualitative differences. As mentioned in HCll, 
these differences arise essentially from the two essential characteristics of our theory, namely: 
(i) each overdense region dominated by gravity collapses under its own, density- dependent 
dynamical time, (ii) there is no particular threshold density/scale for star formation; instead, 
any overdense region produced by initial density fluctuations is susceptible to collapse, no 
matter its degree of internal turbulence, if dominated by gravity. 

We explore the dependence of the SFR upon the clump properties, namely the size, 
the level of turbulence (through the virial parameter) and the magnetic fleld. We show 
that the exact value of the SFR depends signiflcantly upon these clump parameters, which, 
unfortunately, are ill deflned quantities, hampering precise theoretical estimates of the star 
formation rate in a clump/cloud in particular for magnetized clouds. We show that, when 
time-dependence is taken into account in the formation of collapsing dense cores, turbulence 
globally favors star formation, in contrast to the result obtained with the time-independent 
theory, in agreement with the results of Padoan & Nordlund (2011). We also show that, 
when the clump is unbound, i.e. amr > 1, the star formation rate decreases drastically, as 
the dominant contribution of kinetic (thermal+turbulent) over gravitational energy prevents 
gravitational collapse of the overdense structures, inhibiting star formation. The magnetic 
fleld also reduces the SFR when it is strong enough. Its exact influence depends on the 
magnetic-density relation and also on its influence on the density PDF none of these prop- 
erties being known with great accuracy. 

Our calculations show that star formation rate correlates with both the gas surface 
density of the surrounding region and the size of the clump and occurs dominantly in dense 
and/or massive - thus bound - clumps. Such large clumps and dense regions, however, 
represent only a modest fraction (<20%) of typical clouds. In contrast, for the smallest 
clumps, the SFR decreases drastically below a typical density S^, ^ IOOMqPc"^, i.e. a 
volume density n ~ 2500 {Rc/^ pc)~^§, i.e. a visual extinction Ay ~ 7, in excellent agreement 
with observational determinations in Milky Way molecular clouds. For clumps below about 
this density, star formation is still taking place, although at a much lower rate, only in very 
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large (> 10 pc) clumps, which are very rare. This stresses the need to explore large fields 
of view to accurately determine the SFR in low-density regions. We stress that this abrupt 
change in the SFR does not stem from an arbitrarily defined threshold for star formation, 
but naturally arises from the theory and refiects the fact that pc-size low density clumps 
barely form stars essentially because of their too large characteristic Jeans length compared 
with their size and their too low level of turbulence. 

Interestingly enough, the aforementioned density corresponds to the density above which 
the cold gas in the ISM becomes essentially composed of molecular hydrogen, H2. Indeed, 
simulations of compressible turbulence coupled with hydrogen chemistry (Glover & MacLow 
2007, see also Krumholz et al. 2009) suggest that (i) H2 forms much more efficiently in 
(compressible) turbulent gas than in quiescent gas owing to the much shorter formation 
timescale and that (ii) above n ~ 300^^, the H2 formation rate, which scales as the square of 
the number density, n^, becomes larger than the photodissociation rate, allowing the efficient 
in-situ formation of molecular hydrogen. The gas becomes dominantly (resp. entirely) 
molecular above about 5000^ (resp. 10^5) ^"^^ remains "trapped" and shielded from external 
UV radiation in the collapsing structures. Our calculations thus support Glover & MacLow's 
(2007) suggestion that dense (bound or unbound) star-forming regions drive H2 formation, 
and not the opposite. Phrased differently, the formation of molecular hydrogen is not a 
necessary condition for star formation but instead is a consequence of efficient star formation. 

Finally, we derive a fully time-dependent calculation of the core mass spectrum by con- 
sidering a simple time-dependent clump model. That is, instead of assuming fixed clump 
conditions (density, mass, size), we consider an evolving, accreting clump assuming a con- 
stant accretion rate. In order to explore the parameter space, we have conducted calculations 
for various accretion rates, typical of the observed determinations, assuming that the clump 
parameters obey the standard Larson scaling relations. The resulting core mass spectrum 
remains in very good agreement with the observationally derived distribution and repro- 
duces various behaviors observed in simulations of converging fiows. We also investigate 
the core mass spectrum that results from the contribution of a clump distribution such that 
dNc/dMc ex M~^'^. We show that the resulting core mass spectrum is dominated by the 
mass spectrum of the biggest clumps because they contain most of the mass. 

Our SFR model provides the general framework that includes the simplified model 
suggested by Lada et al. (2012), which relies on the fraction of dense (molecular) gas in the 
cloud and thus critically depends on a threshold density. As mentioned above (and shown 
in Fig. [7]), our theory naturally predicts a strong, roughly quadratic correlation between 
the SFR and the gas density, predicted to become dominantly molecular above ~ 5000^1. As 
shown in §4.3. H such a correlation immediately implies a similar correlation with the mass of 
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the clump/cloud itself, i.e. with the mass of molecular hydrogen above the aforementioned 
value. The model suggested by Lada et al. (2012) is thus a direct consequence of our general 
theory of the IMF+SFR. As mentioned above, however, star formation is still predicted 
to take place, although at a much lower rate, in large, low-density clouds, which will be 
composed essentially of atomic hydrogen. 

However, in contrast to what has often been claimed in the literature, we show that 
there is no universal value of the star formation rate, and that this latter does not simply 
correlate linearly with gas density. Indeed, as mentioned above, the SFR strongly depends 
not only on the gas density but also on the clump mass/size, which leads to a large scatter 
in SFR values. 



6.2. A paradigm for star formation and for the theory of the IMF 

The analytical theory described in our previous and present papers, which correctly 
reproduces various observational constraints, as discussed in these papers, suggests the fol- 
lowing paradigm for star formation and for the resulting CMF/IMF: 

1) compressible (shock dominated) large-scale turbulence in the cloud, due to various 
possible mechanisms such as accretion, converging cold and warm flows or star formation 
itself, generates a field of density fluctuations at all scales in the cloud. The (nearly log- 
normal) PDF of these fluctuations is entirely determined by the characteristic (universal) 
log-density power spectrum index of turbulence. This PDF leads to overdense regions which 
correspond to the observed clump mass spectrum (see paper I). Observations suggest that 
the clumps have a filamentary structure. As they accrete mass and dissipate kinetic energy, 
the densest clumps become gravitationally unstable above a typical density Nh2 ~ 10^^ cm~^ 
(i.e. n ~ 10^ (-Rc/1 pc)~^§ for a spherical clump) for a temperature T ~ 10 K, triggering the 
fragmentation into prestellar cores. Although probably quantitatively affecting our results 
by numerical factors of order unity, this filamentary nature of the clumps does not modify 
the general framework of our theory, which indeed predicts a very narrow CMF below about 
this density for pc-size clumps, yielding a negligible number of prestellar cores. 

2) fragmentation then introduces a scale- dependence (thus a density-dependence) in the 
processus, set up the virial criterion for gravitational collapse (see paper I). As demonstrated 
in paper I and in Chabrier & Hennebelle (2011), we suggest that turbulence plays a crucial 
role in setting the massive initial mass reservoir distribution (thus the Salpeter slope of the 
IMF), progenitors of massive stars, above about the mean thermal Jeans mass, at the early 
stages of star formation. Not by providing a pressure support, in a static sense, but by 
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dispersing the gas within these structures, which would otherwise have collapsed, until they 
reach a mass that we identify as their turbulent Jeans mass (see eq. [28]) . 

3) below about the typical density mentioned above, our theory predicts a drastically 
decreasing CMF and thus SFR, for the typical size of the dominant clump population, essen- 
tially because the characteristic Jeans length of the clump becomes comparable to or larger 
than its size, stabilizing the clump against fragmentation. For a filament, the same density 
corresponds to the critical mass per unit length below which the clump is gravitationally 
stable, yielding the same conclusion. We predict, however, that star formation can still occur 
below this density in very large ( > 10)pc size) clumps, large enough to significantly exceed 
their typical Jeans length and generate large enough turbulence levels. This seems to be sup- 
ported by the observed population of class II cores in low-density, large clumps (Gutermuth 
et al. 2011) and stresses the need to observe very large areas at low-density in order to get 
statistically significant core detections. Such clumps, however, are very seldom, as inferred 
from the steeply decreasing clump size spectrum, and the probability to find pretellar cores 
in low-density environments drops drastically, making their detection very difficult. 

4) star formation is thus a continuous process and can occur, statistically speaking, in 
any density environment, as any turbulence- induced overdense region can collapse, if domi- 
nated by gravity, and produce eventually a prestellar core. There is thus no real "threshold" 
for star formation. This leads to a direct correspondence between the star formation rate 
and the gas (atomic -|- molecular) density, as observed and indeed predicted by the theory, 
with star formation occurring most actively in the densest regions of molecular clouds, which 
entails only a small ( < 20%) fraction of their mass. However, as mentioned above, star for- 
mation is basically choked off below n ~ 1000^ for the typical (~pc) size of most clumps. 
Above this density, the H2 formation rate, which increases as the square of the gas density, 
dominates the photo-dissociation rate and the gas quickly becomes dominantly and entirely 
molecular. This drastically decreases the photoelectric heating efficiency, causing in turn a 
drop in dust and gas temperature (Tielens & HoUenbach 1985, Glover & McLow 2007). Star 
formation thus promotes II2 formation, and not the opposite. 

We stress, however, that, from the general point of view, star formation strongly depends 
on the clump characteristic properties (mass, size); therefore, there is no "universal" relation 
between star formation rate and gas density, but instead large variations, depending on the 
clump's environment. 

This paradigm, which relies on the results presented in the present and former papers, 
shows that a star formation rate determined at the early stages by turbulence-induced fluc- 
tuations provides quite a consistent picture of star formation in Milky Way molecular clouds, 
with star formation occurring dominantly in the densest regions of the cloud. This gravotur- 
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bulent picture of star formation correctly predicts the observed CMF over the entire mass 
range from brown dwarfs to massive stars and naturally leads to a star formation rate vs 
clump masses/sizes and gas density correlations in very good agreement with observational 
determinations. We stress, however, that ill-defined quantities such as the size of the clump 
or the dynamical efficiency of core-to-star mass conversion, which may well vary from clump 
to clump, make the determination of the SFR uncertain by a factor of a few. 

The research leading to these results has received funding from the European Research 
Council under the European Community's Seventh Framework Programme (FP7/2007-2013 
Grant Agreement no. 247060. This work was initiated while G.C. was a visitor of the Max 
Planck Institute for Astrophysics, where he benefited from numerous discussions with various 
colleagues. 
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